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1.  INTRODUCTION 


The  U.S.  Air  Force  has  sponsored  several  programs  in  pursuit  of  the 
Automatic  Technical  Control  (ATEC)  of  telecommunications  networks.  Systems 
have  been  developed  which  perform  extensive  "out-of- service"  channel  analysis 
and  which  have  a  limited  capability  for  "in-service”  monitoring.  "Out-of- 
service"  refers  to  channels  not  in  use  at  the  time  of  testing.  "In-service" 
refers  to  non-intrusive  probing  of  a  channel  while  it  is  transferring 
information.  When  a  channel  has  degraded  to  such  a  degree  that  it  will  no 
longer  transfer  data  with  sufficient  accuracy  to  maintain  effective 
communications,  it  has  "failed."  When  a  channel  fails,  it  is  necessary  to 
re-route  the  data  through  a  back-up  channel  and  subject  the  failed  channel  to 
out-of- service  testing.  The  development  of  in-service  channel  monitoring 
techniques  is  desirable  since  It  would  minimize  the  number  of  reserve  channels 
required  and  would  provide  advanced  warning  of  channel  failures  before 
acceptable  service  is  disrupted.  The  implementation  of  such  an  in-service 
monitoring  capability  would  probably  provide  more  reliable  service,  while 
reducing  overall  system  maintenance  costs,  by  scheduling  down  time  on  falling 
channels  for  testing  during  off-peak  hours. 

The  objective  of  this  effort  was  to  theoretically  and  experimentally 
investigate  the  use  of  modem  signature  analysis  techniques  for  in-service 
channel -quality  monitoring  of  voice  bandwidth  channels  used  for  digital 
transmission.  "Modem  signatures”  is  the  term  applied  to  any  characteristic 
derived  from  the  signal  source  by  measurement  or  transformation  (FFT,  for 
example)  for  use  by  a  generalized  algorithm.  Monitoring  refers  to  the 
determination  of  the  transmission  channel  condition  and,  ideally,  to  the 
identification  of  the  type  and  degree  of  impairments  affecting  the  channel,  to 
aid  the  operator  in  effectively  managing  available  resources  (channels)  and 
maximizing  up- time. 
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2.  BACKGROUND 


Recent  developments  In  microprocessor  and  VLSI  technology  have  reduced 
modem  and  computer  costs  and  have  resulted  in  an  increased  utilization  of 
telecommunications  channels  for  data  communication  purposes.  In  response  to 
the  increasing  demand  for  enhanced  in-service  monitoring  capabilities,  the  Air 
Force  Coommlcations  Service  (AFCS),  in  1977.  performed  an  in-house 
investigation  which  showed  that  a  "large  variety"  of  in-service  data  modem 
samples  could  be  distinguished  from  one  another.  After  the  installation  of 
the  Wire-Line  simulator  at  the  Digital  Communications  Experimental  Facility 
(DICEF),  RADC ,  it  became  possible  to  simulate  degraded  channel  conditions  in  a 
controlled  fashion.  With  this  capability,  a  four-modem  database  was 
generated.  The  modems  were  each  processed  through  the  wire-line  simulator 
where  varying  degrees  of  Impairments  were  added.  The  modems  of  the  original 
database  are  listed  in  Table  2-1.  Each  modem  sample  cut  was  approximately  1 
1/2  seconds  long  and  digitized  at  12800  samples  per  second.  The  type  and 
degree  of  impairments  collected  are  shown  in  Table  2-2. 

The  features  used  in  the  experiment  were  functionally  equivalent  to 
calculating  the  autocorrelation  of  the  signal  after  it  had  been  quantized  to  1 
bit.  This  is  shown  graphically  in  Figure  2-1.  Generally,  the  lowest  lags  of 
the  autocorrelation  provided  the  best  discrimination  between  modems. 

It  was  shown  that  the  four-modem  database  could  be  separated  with 
complete  reliability  in  the  presence  of  noise,  harmonic  distortion,  and  phase 
jitter.  No  attempt,  however,  was  made  to  characterize  the  channel  condition. 
The  problem  of  Identifying  signal  featires  and  classification  algorithms  to 
perform  channel  degradation  detection  and  classification  were  the  goals  of  the 
Modem  Signature  Analysis  (MSA)  effort.  The  remainder  of  this  report  discusses 
the  results  of  this  program. 
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3.  DATABASE  OVERVIEW 


Under  the  present  effort,  it  was  agreed  to  expand  the  original  four-modem 
database.  As  more  modems  were  added  to  the  database  it  became  more  difficult 
to  classify  a  modem  type.  Further,  the  problem  of  deducing  the  channel 
quality  had  more  cases  of  "normal"  operating  parameters  to  consider.  It  was, 
however,  decided  that  the  original  four  modems  did  not  provide  a  sufficient 
basis  for  an  investigation  of  this  scope. 

The  seven  additional  modems  collected  under  this  effort,  together  with 
the  original  four ,  are  listed  in  Table  3-1.  All  of  the  later  collections  were 
at  a  sampling  rate  of  17.2  KHz  in  contrast  to  the  earlier  collection  which 
used  12.8  KHz.  All  the  signals  of  the  first  collection  included  cuts  impaired 
with  varying  amounts  of  Gaussian  noise,  harmonic  distortion,  or  phase  Jitter 
at  10,  60  and  120  Hz.  In  addition,  each  signal  included  some  cuts  impaired 
with  a  combination  of  Gaussian  noise  and  phase  jitter  at  60  Hz. 

Under  this  effort,  there  was  concern  that  the  wireline  simulator's 
capability  to  inject  harmonic  distortion  was  not  functioning  correctly. 
Because  of  this,  no  cases  of  harmonic  distortion  were  collected  under  this 
effort.  All  signals  were,  however,  collected  while  impaired  with  Gaussian 
noise  and/or  phase  Jitter.  The  range  of  impairments  collected  for  each  modem 
of  the  complete  database  is  shown  in  Table  3-2. 

The  complete  database  contained  a  total  of  311  modem  samples.  Each  modem 
type  was  supplied  on  a  raw  formatted  tape  accompanied  by  a  log  sheet 
describing  the  30  or  so  cuts  on  that  particular  tape.  The  cross-referencing 
between  the  logsheets  and  data  tapes  soon  became  unmanageable  as  the  database 
grew  larger.  It  was  decided  to  eliminate  the  need  for  the  log  sheet 
referencing  by  encoding  the  modem  type  and  collection  parameters  into  the  name 
of  the  file  which  stored  each  cut  of  data. 
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Software  development  was  performed  on  a  VAX  11/780  and  on  a  POP  11/45. 
The  target  processor  was  a  POP  11/40.  Both  operating  systems  (VAX  VMS  and 
PDP-11  RSX11M)  support  FILES-11  disk  file  structure.  FILES-11  allows  each 
file  to  be  named  by  a  nine-character  name  and  a  three-character  extension. 

The  filename  encoding  scheme  employed  is  shown  in  Table  3-3.  The  first 
four  characters  include  the  generic  code  information,  including  the  number  of 
bits/ second,  the  baud  rate,  and  modulation  type.  The  first  character  is  an 
arbitrary,  unique  symbol  assigned  to  represent  each  Individual  modem  type. 
Table  3-1  includes  the  generic  code  label  for  each  modem. 

Several  benefits  are  immediately  evident  from  this  scheme.  For  example, 
if  all  modems  with  a  baud  rate  of  2400  are  desired,  it  is  simply  necessary  to 
scan  all  the  file  names  for  those  with  a  "3"  in  the  third  position. 
Furthermore,  reviewing  the  success  of  a  processing  run  is  simplified  because 
all  available  groundtruth  is  encoded  into  the  filename  of  the  data  sample. 

In  addition  to  the  generic  and  ID  information  encoded  into  the  filename, 
the  impairments  present  on  the  channel  are  also  encoded.  The  fifth  character 
describes  the  nature  of  the  impairment.  As  shown  in  Table  3-3 •  codes  "0”  and 
"I"  refer  to  the  raw  modem  signal  which  has  not  been  passed  through  the  4A 
line,  and  to  the  modem  signal  which  has  been  passed  through  the  4A  line  with 
no  impairments  added,  respectively.  The  two-character  field  "FI"  identifies 
the  amount  of  Gaussian  noise  or  harmonic  distortion  added  to  the  sample  in 
dB's.  Field  "F2"  indicates  the  amount  of  phase  Jitter  in  degrees  peak-to- 
peak.  This  allows  the  encoding  of  combined  noise  and  phase  Jitter  information 
in  the  filename.  Note  that  whenever  phase  Jitter  occprs  in  combination  with 
Gaussian  noise,  the  frequency  of  the  phase  Jitter  component  is  always  60  Hx. 
This  case  corresponds  to  an  "I*  code  of  "4".  The  other  frequencies  of  phase 
Jitter  each  are  identified  by  a  separate  "I"  code  of  "S"  through  "S";  field 
"F2”  shows  the  degrees  peak-to-peak. 
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Table  3-3 


MODEM  ID  DATA  BASE  FILENAME  CONVENTION 


T 

#  # 

N 

* 

I  Fl  F2  " 

HAV" 

/ 

TYPE 

BITS/SEC. 

BAUD 

s 

MODULATION 

p 

- 

PARADYNE  9600 

1-1200 

1-1200 

F-FSK 

2-2400 

2-1600 

D- PSK/DPSK 

L 

- 

LENKURT  26C 

3-4800 

3-24Q0 

Q-QAM 

4-9600 

4-2667 

P-PAM 

H 

- 

HUGHES  276 

S-16000 

5-4800 

C 

CODEX  LSI  9600 

(I)  IMPAIRMENT  (PI) 

FIELD  1  (F2) 

FIEDD  2 

M 

_ 

PARADYNE  MP 

t  -  NO  4A 

00 

00 

1  -  4A  CLEAN 

00 

00 

A 

- 

CODEX  LSI  48  MODE 

A 

2  -  HARMONIC  DIST. 

+  dB 

00 

3  -  GAUSSIAN  NOISE 

♦  dB 

00 

B 

- 

CODEX  LSI  48  MODE 

B 

4  -  NOISE 

(  PHASE  JITTER  +dB 

P/P® 

5  -  CALIBRATION  SIGNAL  00 

00 

D 

- 

CODEX  LSI  48  MODE 

C 

6  -  JITTER 

0  10  hr 

00 

P/P° 

7  -  JITTER 

0  60  hz 

00 

P/P° 

W 

- 

WECO 

8  -  JITTER 

0  120  bz 

00 

P/P® 

9  -  JITTER 

9  180  bz 

00 

P/P° 

R 

- 

HARRIS 

X 

- 

MD-674 

Table  3-3 
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Code  "5",  described  as  calibration  signals  in  Table  3-3.  includes  a 
variety  of  test  signals.  These  are  typically  test  tones  or  Gaussian  noise 
passed  through  the  simulated  4A  channel. 

3.1  DATA  COLLECTION 

All  of  the  data  collection  was  performed  at  the  RADC/DICEF,  Griffiss  Air 
Force  Base  by  RADC  personnel.  The  original  data  collection  setup  is  shown  in 
Figure  3-1.  The  random  sequence  generator  provides  a  repetitive  random  stream 
into  the  modem  under  test.  The  modem  output  is  fed  to  the  wireline  simulator. 
The  signal  is  digitized  in  the  simulator  and  the  channel  characteristics  and 
programmed  impairments  are  added.  The  resultant  signal  is  D/A  converted  and 
then  leaves  the  simulator.  After  anti-aliasing  filtering  the  signal  is  re¬ 
digitized  at  12.8  KHz  for  data  collection  and  recorded  on  7-track  magtape  by 
the  DICEF  9303  processor.  For  the  purpose  of  this  investigation,  the  data  was 
converted  by  RADC  personnel  to  9-track  magtape  and  supplied  to  PAR  Technology 
Corporation. 

The  original  database  for  the  MSA  effort  contained  4  modem  types: 

•  CDX  LSI-9600  -  QAM 

•  Hughes  HC-276  -  DPSK 

•  Paradyne  LSI-96  -  PAM  (Vestigual  Side  Band) 

•  Lenkurt  26-C  -  FSK 

The  second  series  of  data  collections  were  performed  on  the  setup  shown 
in  Figure  3 .2  Note  the  51.6  KHz  sample  rate  into  the  AP120B.  Data  was 
resampled  by  three  to  reduce  the  sampling  rate  to  17.2  KHz  before  being  passed 
to  the  PDP  11/40  processor  and  saved  to  magtape. 
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Figure  3-2  Second  MSA  Data  Collection  System 


From  the  two  collections  a  total  of  11  modems  were  sampled.  Four  sampled 
at  12.8  KHz  and  seven  sampled  at  17.2  KHz.  In  order  that  features  extracted 
from  the  modem  signals  be  comparable,  it  is  required  that  all  of  the  samples 
have  the  same  sampling  rate.  Therefore,  it  was  determined  that  the  sample 
rate  of  the  database  be  unified.  Further,  in  order  that  minimum  error  be 
introduced  into  the  signals,  it  was  decided  that  the  17.2  KHz  sampling  rate 
would  be  resampled  down  to  approximate  the  lower  sampling  rate,  12.8  KHz. 

A  method  known  in  engineering  as  Sampling  Rate  Decimation,  was  applied  to 
the  17.2  KHz  sampled  selected  Modem  types  in  the  MSA  database  to  alleviate 
this  sampling  rate  discrepancy.  The  implementation  of  this  sampling  reduction 

algorithm  is  fairly  straight  forward  and  is  described  elsewhere1.  A  summary 
of  the  method  is  given  here. 

The  two  sampling  rates  contained  in  the  MSA  database  were  12.8  KHz  and 
17.2  KHz.  Since  12.8  KHz  is  approximately  three-quarters  of  17.2  KHz,  the 
sampling  rate  reduction  technique  to  be  applied  is  straightforward.  The 
approach  accomplishes  the  reduction  of  the  sampling  rate  from  17.2  KHz  to  12.9 
KHz  within  0.81  of  the  desired  sampling  rate  12.8  KHz.  The  procedure  is 
indicated  below: 

1.  Initial  Data  is  obtained  at  a  sampling  rate  of  17.2  KHz. 

2.  Pad  each  sample  in  the  Raw  Input  Data  File  with  2  zeros 

3.  Apply  a  Low  pass  filter  with  a  6  KHz  cutoff  and  a  gain  of  3.0 

4.  Resample  by  a  factor  of  4 


1  Proceeding  of  the  IEEE,  Vol.  69,  No.  3.  March  1981  "Interpolation  and 
Decimation  of  Digital  Signals  —  A  Tutorial  Review."  by  Ronald  E. 
Crochlere,  Lawrence  R.  Rablner. 
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After  execution  of  step  (2),  the  data  stream  will  have  an  effective 
sampling  rate  of  51.6  KHz.  The  6  KHz  low-pass,  linear-phase  FIR  filter  will 
Interpolate  the  padded  zeros  and  act  as  an  anti-aliasing  filter. 

Since  the  effective  sampling  rate  is  currently  51.6  KHz,  resampling  of 
the  data  by  a  factor  of  4  yields  a  sampling  rate  of  12.9  KHz,  which  is  within 
0.81  of  our  desired  sampling  rate,  12.8  KHz.  Following  the  conversion  of  all 
modem  types  originally  collected  at  a  sampling  rate  of  17.2  KHz  to  the  12.9 
KHz  sample  rate,  extensive  tests  were  performed  to  Insure  the  validity  of  the 
resampled  data.  Figure  3-3  shows  a  time  plot  containing  the  first  1024 
samples  of  the  PARADYNE  MP-96  modem  with  AGN  added. 

This  illustration  was  plotted  from  the  original  modem  database  which  was 
generated  at  a  17.2  KHz  sampling  rate.  Figure  3-4  illustrates  a  power 
spectral  density  plot  of  the  same  data.  Figure  3-5  shows  the  same  data  after 
resampling  (sampling  rate  =  12.9  KHz).  The  difference  between  the  plot  shown 
in  Figure  3-5  and  the  plot  shown  in  Figure  3-3  indicates  the  error  Introduced 
by  the  reduction  of  the  sample  rate.  By  visual  Inspection,  it  is  clear  that 
the  plots  are  identical.  Inspection  of  the  Power  Spectral  Densities  of  the 
same  modem  data  at  different  sampling  rates,  Figures  3-4  and  3-6,  we  concluded 
that  the  carrier  frequencies  and  bandwidths  of  both  signals  are  equal,  thus 
showing  the  success  of  the  sampling  rate  reduction. 

The  sampling  reduction  technique  was  applied  to  each  file  in  the  MSA 
database  which  has  been  generated  using  a  sampling  rate  of  17.2  KHz.  The 
resulting  database  currently  exists  with  a  sampling  rate  of  12.9  KHz. 

Figure  3-7  illustrates  the  frequency  response  of  the  6  KHz  low-pass 
filter  used  in  the  resampling  technique.  Figure  3-8  lists  the  coefficients 
which  make  up  the  6  KHz  filter. 


Time  Plot  Paradyne  MP-96  AGN  20  dB  S/N 


DC  Removed,  Hanning 


1620  Wl 


258 

Figure  3-6:  Power  Spectral  Density,  Paradyne  MP-96,  AGN  20  dB  S/N  5  Line  Avg. 

DC  Removed,  Hanning  Weighting,  60  dBi  fcg  Scale 


QUANTIZATION  IN  BITSi  22  GAIN  IN  OBi  *.5*  NUMBEA  OF  TAPBi  64 
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Figure  3-7:  Frequency  Response  for  6kHz  Lowpass  Filter  Used  in  Sample  Reduction 

Algorithm  Implementation  Filter  Length  *  64  Taps 
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3-8:  Top  Weights  of  the  Lowpass  Filter  Used  for  Resampling 


3.2  THE  WIRELINE  SIMULATOR 


The  wireline  simulator  was  used  to  produce  all  of  the  channel  Impairments 
which  were  compiled  into  the  MSA  database.  Since  the  MSA  database  will  be 
studied  to  suggest  possible  algorithms  and  will  also  serve  as  a  test  data  set 
for  any  algorithms  developed,  how  the  simulator  performs  the  corruption  of  the 
channel  will  be  described.  The  description  will  focus  on  the  theory  of 
operation  rather  than  the  implementation  details  and  will  emphasize  the 
particular  parameters  used  during  this  program  rather  than  the  myriad  of 
options.  For  further  information  the  reader  is  referred  to  at  least  three 
detailed  technical  manuals  available  in  the  DICEF  on  this  subject. 

The  wireline  simulator  was  originally  installed  at  DICEF  in  1970.  In 
1976  the  simulator  was  upgraded  with  the  addition  of  a  special  purpose 
programmable  signal  processor  and  a  new  control  minicomputer.  The  description 
can  be  conveniently  separated  into  the  original  simulator  and  the  upgrade.  A 
block  diagram  of  the  original  is  shown  in  Figure  3.9.  It  comprises  a 
programmable  digital  filter,  a  frequency  offset  unit,  an  Impulse  noise 
generator,  and  interfaces  with  a  general  purpose  digital  computer  and  with 
input  and  output  analog  signals. 

The  digital  filter  is  based  on  a  direct  convolutional  algorithm  employing 
a  special  high-speed  arithmetic  unit  with  a  400  nanosecond  multiply/add  time. 
The  frequency  offset  unit  performs  a  digital  single- sideband  modulation 
including  jitter  of  the  offset  frequency.  The  impulse  generator  produces 
rectangular  pulses,  noise  bursts  and  damped  sinusoids  each  with  variable 
parameters.  Automatic  amplitude  scaling  A/D  and  D/A  converters  provide  an 
Input/output  dynamic  range  well  in  excess  of  the  basic  converter  accuracy. 

A  block  diagram  of  the  upgraded  wireline  simulator  is  shown  in  Figure  3- 
10.  In  this  figure  the  original  simulator  is  represented  by  the  large  block 
in  the  left  center.  It  is  noted  that  the  upgrade  was  added  in  a  serial 
fashion  to  the  original.  That  is,  before  digital-to-analog  conversion  at  the 
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Figure  3-9:  Original  Wireline  Simulator 


Figure  3-10:  Upgraded  Wireline  Simulator 


output  port  of  the  original,  the  sampled  signal  was  Interrupted  and  sent  to 
the  upgrade  for  the  addition  of  more  parameters. 

The  first  parameter  of  interest  is  additive  Gaussian  noise.  An  external 
General  Radio  Type  1390-A  Random  Noise  Generator  was  applied  to  the  noise 
input  shown  in  Figure  3-9.  This  source  was  applied  before  the  input  A/D 
converter  and  as  a  result  was  weighted  by  the  remainder  of  the  simulation. 

The  next  major  component  applicable  to  this  program  is  the  programmable 
digital  filter  used  to  simulate  the  attenuation  and  delay  distortion 
characteristics  of  the  wireline  channel.  The  impulse  response  is  represented 
by  256  samples.  Thus,  each  output  sample  is  generated  by  256  multiply/ add 
operations  in  the  convolutional  filter.  The  filter  performs  a  multiply/ add 
operation  in  400  nanoseconds.  This  results  in  a  sampled  rate  of  1/(400  x  10“^ 
x  256)  a  9765.50  samples/  sec.  The  filter  employed  a  typical  convolution 
algorithm  where  the  nth  sample  of  the  output  function  g(n),  is  given  by 

N 

g(n)  *  £  h  (k)  f  (n  -  k)  (1) 

k*0 

where  h(k)  is  the  sampled  impulse  response  of  the  filter  and  f(k)  is  the  input 
function.  Implementation  of  Eq.  (1)  requires  two  circulating  memories,  one 
which  stores  h(k)  and  one  which  stores  f(n  -  k) ,  and  a  multiplication/ addition 
unit  which  generates  g(n).  Note  that  as  n  advances,  a  particular  sample  of 
f(k  -  n)  slides  relative  to  h(k).  This  is  achieved  in  practice  by  making  the 
memory  which  stores  f(n  -  k)  one  word  longer  than  the  memory  which  stores 
h(k).  Then,  when  the  h(k)  memory  completes  a  full  cycle  and  a  new  value  of 
g(n)  has  been  generated,  the  f(n  -  k)  memory  lacks  one  word  of  completing  a 
full  cycle.  Thus  the  data  in  the  f(n  -k)  memory  is  shifted  back  by  one 
position  relative  to  the  data  in  the  h(k)  memory  every  cycle.  The  vacant 
position  at  the  beginning  of  the  f(n  -  k)  memory  is  filled  by  the  new  input 
sample. 
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Figure  3-11  shows  a  block  diagram  of  the  programmable  digital  filter  used 
in  the  simulator.  The  hCk)  memory  contains  256  samples  of  the  filter  impulse 
response  and  the  f(n  -  k)  memory  contains  257  samples  of  the  input  function. 
The  particular  impulse  response  required  for  the  filter  is  loaded  into  the 
h(k)  memory  from  paper  tape,  through  a  paper  tape  reader  and  control  unit. 

Voice  channel  communications  often  go  through  carrier  systems  which 
employ  single  sideband  modulation  techniques.  Usually,  the  transmitter  and 
receiver  carriers  are  not  exactly  locked.  In  particular,  one  carrier  will 
drift  in  frequency  and  phase  relative  to  the  other  one,  producing  a  variable 
spectrun  shift  in  the  data  signal.  The  frequency  shift  unit  is  Intended  to 
simulate  this  effect.  It  is  capable  of  introducing  a  frequency  offset  in  the 
spectrun  of  the  data  signal.  The  frequency  offset  can  be  modulated  by  either 
a  square  or  sinusoidal  modulation  waveform  where  the  frequency  and  deviation 
of  the  modulation  can  be  controlled.  This  unit  was  used  to  impart  phase 
jitter  to  the  MSA  database. 

The  frequency  shift  unit,  as  Implemented,  performs  a  single  sideband 
modulation  directly  at  audio  frequency.  The  mathematical  block  diagram  of  the 
unit  is  shown  in  Figure  3-12  and  a  functional  diagram  is  shown  in  Figure  3-13* 
The  mathematics  of  the  signal  processing  are  outlined  below. 

The  input  to  the  unit,  f(t),  is  to  be  shifted  in  frequency  byOJ^.  Since 
the  duration  of  f(t)  is  always  finite,  it  can  be  represented  over  its  entire 
duration  by  a  Fourier  series 

f(t)  -  £  [a  c«.  s  n^t  +  sin  n(j£  t]  (2) 

n 


The  signal  f(t)  is  first  passed  through  a  90°  phase  shift  device  and  the 
resulting  quadrature  signal  is  fq(t),  where 

f  q  ( t )  -  £  [bn  cos  n(j^t  -  an  sin  n(^t]  (3) 

n 
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Figure  3-11:  Digital  Filter  Block  Diagram 


Digital  Frequency  Shift  Unit 


Finally,  the  single  sideband  modulation  is  performed  as  the  difference  of  two 
products : 


fg(t)  -  fq(t)  sinC0ot  -  f(t)cos(|)ot 

-E  [bn(:os  nC*V  sil**V  +  8itl  nW0t  C°8<*V} 

li  9 

+  an  [cos  nWot  cosCOmt  -  sin  n^t  sin^t] 

~L[\  cos(nW0  “^m>c  +  bn  sin(na)o  -Um)l] 


Thus,  fs(t)  is  the  same  signal  as  f(t)  with  the  spectrum  shifted  by  an 

increment  o)_. 

m 

The  next  impairment  function  to  be  considered  is  the  harmonic  distortion 
generator.  This  device  presented  a  unique  design  problem  since  for  high 
distortion  levels  it  is  possible  to  produce  spectral  widening  with  its 
associated  aliasing.  The  inputs  to  the  distortion  generation  elements  are 
band  limited  so  that  the  distorted  output  signal  will  not  contain  significant 
aliasing.  This  is  accomplished  as  shown  in  Figure  3-14.  All  low-pass  filters 
are  compensated  for  zero  delay  and  zero  delay  distortion. 

The  input  signal  which  is  sampled  at  9765.63  Hz,  is  first  band  limited  to 
2500  Hz  and  then  squared.  The  output  of  the  squarer  will  contain  second 
harmonic  components  which  have  a  maximum  frequency  at  5000  Hz.  Thus,  no 
lnband  aliasing  is  produced.  That  output  will  in  turn  be  band  limited  to  2500 
Hz  and  again  squared  to  produce  components  with  4th  harmonic  components.  The 
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Distortion  Generator 


second  and  the  filtered  4th  harmonic  components  are  multiplied  to  produce  6th 
harmonic  components.  The  odd  harmonic  components  are  generated  by  mixing  even 
harmonic  components  with  the  fundamental.  Examination  of  the  figure  will 
reveal  that  no  signal  component  used  will  have  a  spectrin  wider  than  half  the 
sampling  rate.  This  will  avoid  all  aliasing  problems.  It  should  be  noted 
that  none  of  these  components  is  pure,  but  the  appropriate  harmonic  energy  is 
present  and  by  combination  of  the  components,  it  is  possible  to  generate 
arbitrary  harmonic  distortion.  This  capability  is  provided  by  the  weighting 
network  which  implements  the  coefficients  of  a  Taylor  series  expansion  of  the 
distortion  characteristic  to  produce  the  desired  harmonic  energy.  The 

operator  enters  harmonic  amplitude  and  the  PDP-8  computes  the  coefficients. 
When  the  harmonic  distortion  signal  has  been  generated,  the  result  is  then 
added  to  the  original,  undistorted,  signal. 

The  harmonic  distortion  process  is  performed  by  combining  the  terms  in  a 
Taylor  series  expansion. 


Distortion  Component  *  aQ  +  a-,  S(t)  +  a2S2(t)  +  a3s3(t) 

+  a^s^Ct)  ♦  8gS5(t)  a6s6(t) 


where  S(t)  is  the  input  signal  to  the  distortion  generator.  The  6th  order 
term  in  the  expansion  also  contains  4th  and  2nd  harmonic  energy  in  addition  to 
the  6th  harmonic  energy.  The  PDP-8  OPCOM  program  is  designed  to  operate  on 
operator  inputs  specifying  harmonic  levels  to  create  the  set  of  coefficients 
for  the  Taylor  series  which  cancel  undesired  harmonics.  Note,  this  procedure 
can  function  only  when  signal  level  is  known  and  S(t)  is  a  sine  wave.  It  is 
possible  to  generate  pure  6th  harmonic  by  requesting  the  OPCGM  program  to 
generate  only  a  6th  harmonic  component.  The  3rd  harmonic  was  utilized  for  the 
MSA  data  collection. 
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4.  A  HIERARCHICAL  APPROACH  TO  PERFORMANCE  MONITORING  OF  VOICE  FREQUENCY  CHANNELS 

It  is  desired  to  develop  simple  and  reliable  techniques  for  performance 
monitoring  of  voice  frequency  communications  channels.  The  problem  is 
considerably  simplified  if  the  performance  monitoring  unit  can  observe  modem 
outputs.  In  particular,  a  lumber  of  existing  techniques  are  available  in  this 
case  for  predicting  bit  error  probability,  which  is  often  the  single  important 
measure  of  channel  quality.  Unfortunately,  in  the  present  application  it 
cannot  be  assumed  that  the  modem  outputs  are  available.  This  severely 
complicates  the  problem  and  requires  a  careful  hierarchical  approach  if  the 
design  objectives  are  to  be  achieved.  In  this  section,  such  an  approach  is 
outlined . 

The  decision  hierarchy  is  illustrated  in  Figure  4-1.  On  the  basis  of 
observations  taken  on  the  voice  frequency  (VF)  line,  and  without  knowledge  of 
the  modem  type,  a  preliminary  and  obviously  gross  assessment  must  be  made 
concerning  whether  or  not  the  channel  has  suffered  a  catastrophic  failure. 
This  might  be  the  result  of  total  loss  of  signal,  or  even  the  SNR  ratio 
falling  to  an  unacceptably  low  value.  Here  in  the  first  step  (Phase  1)  of  the 
hierarchy,  some  tests  dealing  with  estimates  of  the  channel  SNR,  or  possibly 
tests  using  spectral  analysis  techniques  may  be  considered. 

Assuming  that  there  is  no  indication  of  gross  channel  failure,  the  next 
phase  (2)  is  the  task  of  specifically  identifying  the  modem  type  in  service  on 
the  VF  communication  channel.  This  level  of  identification  is  necessary  in 
order  to  be  able  to  predict  what  is  the  "normal"  channel  activity  for  this 
modulation  type.  The  classifier  need  not  be  completely  effective.  The 
classifier  may  fail  to  identify  the  correct  modulation  format  with  a 
sufficient  degree  of  certainty.  In  this  case  the  classifier  may  be  directed 
to  repeat  the  classification  attempt  on  successive  channel  samples  a  fixed 
number  of  times.  If  a  confident  decision  has  not  been  achieved  in  this 
period,  it  can  be  declared  that  a  degradation  in  the  modem  signal  structure 
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has  occurred,  rendering  the  signal  unrecognizable,  or  that  the  actual  modem 
type  is  not  a  member  of  the  established  set  which  the  classifier  is  prepared 
to  recognize. 

Assuming  the  confident  classification  of  the  modulation  format  in  Phase 
2,  it  is  then  required  to  provide  some  quantitative  evaluation  of  the  amount 
and  type  of  degradation  affecting  the  channel.  Phase  3  of  the  hierarchy  is 
the  quantitative  evaluation  of  degradation  and  is  by  far  the  most  difficult  of 
the  hierarchical  tasks  to  be  performed.  It  may  require  replication  of  typical 
modes  receive  functions  such  as  bit  synchronization,  phase  tracking,  data 
demodulation,  and  in  some  cases  adaptive  equalization.  It  should  be  noted 
that  the  actual  implementation  of  the  modem  receive  functions  may  differ  from 
that  of  any  specific  modem  types.  That  is,  if  a  correct  identification  of  the 
modulation  format  has  been  made  in  Phase  2,  then  it  is  possible  in  principle 
to  construct  (in  software)  a  generic  modem  which  is  functionally  equivalent  to 
the  actual  modem  on  the  VF  line  under  test.  While  in  principle  this 
replication  of  modem  receive  functions  is  possible,  it  is  potentially  complex 
in  terms  of  implementation.  Nevertheless,  we  proceeded  under  the  assumption 
that  any  implementation  can  be  achieved,  so  that  the  direction  of  the 

Investigation  would  not  be  biased  toward  the  readily  realizable  solution. 


5.  ALGORITHM  INVESTIGATION 


Throughout  the  course  of  this  effort,  a  variety  of  algorithms  were 
investigated  for  all  phases  of  the  hierarchical  approach  described  in  Section 
4.  It  is  frequent  that  algorithms  overlap  the  boundaries  between  hierarchical 
tasks,  being  able  to  contribute  to  more  than  one  phase  of  the  problem.  For 
this  reason,  the  algorithms  investigated  will  be  presented  here  without  strict 
observance  in  identifying  the  precise  hierarchy.  Rather,  the  choice  of 
algorithms  to  be  demonstrated,  and  the  structure  of  the  demonstration  system 
will  be  the  subject  of  Section  6. 

5.1  MODEM  SPECTRAL  PROFILE  ANALYSIS 

The  goal  of  the  spectral  analysis  Investigation  was  to  evaluate  the  use 
of  signal  power  spectra  for  use  as  a  feature  space  from  which  decisions 
concerning  the  channel  traffic  type  and/or  channel  quality  could  more  readily 
be  performed.  One  member  of  each  modem  group  was  selected.  The  spectral 
estimates  of  these  signals  are  shown  in  Figures  5-1  through  5-11.  Each  plot 
represents  the  average  of  ten  1024>point  power  spectra.  Hanning  weighting  has 
been  applied  and  the  plot  is  clipped  off  40  dB  down  from  the  spectrum  maximun. 
The  horizontal  divisions  each  mark  10dB.  The  frequency  range  is  from  0  to  6400 
Hz,  or  256  Hz  per  vertical  tic  mark.  Any  DC  component  has  been  removed. 

The  first  two  spectra  illustrate  the  FSK  modulation  type.  The  familiar 
carrier  frequencies  are  easily  detectable.  The  three  DPSK,  1200  baud  modems 
are  shown  in  Figures  5-3  through  5-5.  The  lower  edge  of  the  pass  band  region 
for  each  sample  contains  spurious  channel  data,  such  as  noise.  The  roll-off 
of  the  lower-edge  for  the  three  modems  is  not  identical,  which  is  required  if 
any  feature  measurements  are  to  be  made  in  this  region  for  the  purposes  of 
generic  type  identification.  A  potentially  valuable  feature  is  shown  in  the 
PAM,  4800  baud  signal  spectrimi  (PARADYNE  LSI-96)  shown  in  Figure  5-10.  Here 
the  2853  Hz  carrier  is  apparent  above  the  other  pass  band  components. 
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LENKURT  26-C  Signal  Spectrum  Average  of  10  Lines,  Hanning 
Weighted  on  a  40  dB  Log  Scale 


Figure  5-3:  HUGHES  HC-276  Signal  Spectrum  Average  of  10  Lines,  Hanning 

Weighted  on  a  40  dB  Log  Scale 


96  Signal  Spectrum  Average  of  10  Lines,  Hanning 
Weighted  on  a  40  dB  Log  Scale 
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Figure  5-10:  PARADYNE  LSI-96  Signal  Spectrum  Average  of  10. Lines,  Hanning 

Weighted  on  a  40  dB  Log  Scale 


Another  technique  used  to  reveal  signal  properties  was  to  apply  full-wave 
rectifioatlon  prior  to  the  speetriai  computation.  This  process  was  shown  to  be 
very  effective  in  enhancing  the  spectral  component  at  the  baud  rate.  Figures 
5-12  through  5-15  show  the  power  speotral  average  for  the  FSK  1200  baud,  DPSK 
1200  baud,  DPSK  1600  baud  and  QAM  1600  baud  signals,  respectively.  No  baud 
component  enhancement  occurred  for  the  higher  speed  modems. 

5.2  AUTO-CORRELATION  ANALYSIS 

As  a  prerequisite  to  the  determination  of  the  type  and  degree  of 
impairments  affecting  the  channel  under  study,  an  identification  of  the 
channel  traffic  to  a  level  of  generic  type  must  be  performed.  Without  such 
information,  it  is  not  known  what  signal  parameters  represent  a  "normal”  low- 

error  operating  mode. 

Previous  attempts  at  classification  of  the  modems  used  an  equivalent  to 
the  autocorrelation  function  generated  from  a  1-bit  quantized  version  of  the 
modem  signal.  Based  on  the  previous  success  of  these  features,  the 
autocorrelation  function  was  used  as  the  feature  space  in  which  the  generic 
modem  type  separation  would  be  investigated.  The  classification  problem  has 
become  greatly  complicated  since  the  work  performed  on  the  original  four  modem 
database.  The  11  modems  studied  represent  seven  unique  generic  signaling 
strategies.  The  goal  of  this  investigation  is  to  design  a  practical 
classifier  for  the  seven  generic  classes  to  determine  the  modulation  type, 
baud  rate,  and  carrier  location  information  required  in  the  forthcoming 
channel  impairment  characterization. 

The  training  set  for  the  classifier  corresponded  to  the  first  30  lags  of 
the  autocorrelation  function  given  by 

M-m-1 

=  I  A(m+q)  *  A(q) 

q  =  ♦ 
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Spectral  Average  of  Rectified  Signal  on 
i  40  dB  Log  Scale 


where: 


m  is  the  desired  lag 

M  is  the  number  of  data  points  in  the  correlation  window 
A  is  the  input  modem  signal  waveform 


Clearly,  lag  (0)  corresponds  to  the  input  signal  squared,  or  the  total 
power  for  the  zero-mean  input  signal .  In  order  to  minimize  the  effects  of 
collection  gain  variations,  lag  (0)  was  always  used  to  normalize  the  feature 
vector  (of  the  first  30  lags). 

After  much  preliminary  investigation,  it  was  decided  to  narrow  the  scope 
of  the  investigation  to  two  autocorrelation  window  lengths,  3,000  samples  and 
8,000  samples.  Both  the  un-quantlzed  signal  and  the  1-bit  quantized  signal 
were  considered;  where  the  1-bit  quantization  is  defined  as: 

IF  A  (i)  >  0.  then  A  ( i)  =  1 ;  else  A  (i)  =  SS 

It  was  concluded  that  both  the  1-bit  quantized  signal  waveform  autocorrelation 
and  the  full  12-bit  quantized  version  showed  merit.  These  four  cases  will  be 
explored  in  the  remainder  of  this  section. 

Figures  5-16  through  5-26  show  a  plot  of  the  autocorrelation  function  of 
each  modem  for  the  case  of  an  8,000  sample  correlation  window  of  1-bit 
quantized  signal  data.  Also  contained  in  these  figures  is  the  signal  spectrum 
estimate  obtained  from  performing  an  FFT  on  the  autocorrelation  function. 
During  this  investigation,  only  the  autocorrelation  functions  were  used  as 
features.  The  signal  spectrum  is  provided  for  inspection  purposes  only. 
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o-Correlation  Function  for  the  HUGHES 
276  Modem 
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Correlation  Function  for  the  CODEX 
B  Modem 


Figure  5-21:  First  64  Lags  of  the  Auto-Correlation  Function  for  the  CODEX 

LSI -48  Mode  A  Modem 


5-23:  First  64  Lags  of  the  Auto-Correlation  Function  for  the  CODEX 

LSI-96  Modem 


Figure  5-24:  First  64  Lags  of  the  Auto-Correlation  Function  for  the 

PARADYNE  LSI-96  Modem 


of  the  Auto-Correlation  Function  for  the 
PARADYNE  MP-96  Modem 


Figure  5-26:  First  64  Lags  of 


5.2.1  OLPARS  Investigation  Of  A  Practical  Generic  Classifier 


The  pattern  recognition  problem  of  the  MODEM  Signature  Analysis  effort  is 
composed  of  feature  extraction,  that  is.  the  definition  of  measurements,  and 
of  pattern  classification.  The  objective  in  selecting  features  for  the  MSA 
data  base  is  to  provide  a  set  of  measurements  which  yield  information  to  aid 
in  discriminating  among  the  various  modem  generic  types,  even  in  the  presence 
of  severe  channel  distortions.  The  features  used  for  discriminating  among 
different  modem  types  in  the  OLPARS  analysis  are  the  lag  amplitudes  of  the 
autocorrelation  function  of  the  raw  time  waveforms  and  the  1-blt  quantized 
waveforms  of  the  MSA  data  base.  The  analysis  of  the  MSA  data  base  was  aided 
by  use  of  the  OLPARS  system  which  helped  perform  both  the  pattern  analysis  and 
pattern  classification  for  the  MSA  effort. 

After  computation  of  the  feature  vectors  from  the  autocorrelation 
functions,  we  were  ready  to  begin  a  pattern  recognition  and  classification 
analysis  on  the  data.  The  OLPARS  system  uses  a  structured  vector  filing 
system  for  data  handling.  A  tree  structure  was  created  fhom  the  MSA  data  base 
which  included  7  classes,  each  class  being  a  particular  modem  generic  type. 
Figure  5-27  illustrates  the  7-class  tree  structure.  Each  class  such  as 
C1200PSK  and  C1600PSK  represent  the  modems  with  modulation  strategies  of  PSK 
at  1200  and  1600  baud,  respectively.  Following  each  class  name  is  the  assigned 
class  symbol  and  the  total  number  of  vectors  for  each  sub-class  are  printed. 

The  modem  naming  convention  seen  throughout  the  OLPARS  analysis  is  that 
presented  in  Table  3-3.  The  right-most  column  of  names  in  Figure  5-27  is  the 
list  of  vector  file  names.  Each  of  these  names  represents  one  modem  type  for 
all  cases  of  particular  impairment.  For  example,  the  first  vector  file  name 
of  the  right-most  column  of  Figure  5-27  is  "LI IF 2".  Table  3-3  decodes  this  as 
"LENKURT  26-C,  1200  bits  per  second,  1200  baud,  FSK,  harmonic  distortion 

added".  This  file  contains  all  cases  of  harmonic  distortion  for  the  LENKURT 
26-6  modem.  Vector  file  name  "L11F3"  decodes  as  "LENKURT  26-C,  1200  bits  per 
second,  1200  baud,  FSK,  Gaussian  noise  added",  and  contains  all  cases  of 
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iyure  S-27 :  Data  Tree  Structure  Illustration  for  the  Generic  Modem  Type  Separability  lnvestiyat ion 


additive  Gaussian  noise  for  the  LENKURT  26-C  modem.  "L11F6"  contains  all 
cases  of  added  phase  jitter  for  the  LENKURT  26-C  modem.  Note  that  for  the 
OLPARS  analysis,  all  cases  of  added  phase  jitter  for  a  particular  modem, 
regardless  of  the  frequency,  were  combined  into  a  single  vector  file  and 
labeled  by  the  impairment  code  "6".  (Table  3-3  identifies  impairment  code  "6" 
as  10  Hz  phase  jitter).  Next  vectors  files  L11F2,  L11F3  and  L11F6  which 
contain  cases  of  harmonic  distortion,  Gaussian  noise  and  phase  jitter, 
respectively,  are  combined  into  "LI IF"  (on  Figure  5-27).  "LI IF"  decodes  as 
"LENKURT  26-6  1200  bits  per  second,  1200  baud,  FSK"  and  contains  all  cases  of 
all  Impairments  for  the  LENKURT  26-C  modem.  Similarly,  "XI IF"  contains  all 
cases  of  all  available  impairments  for  "MD-671*,  1200  bits  per  second,  1200 
baud,  FSK".  "LI IF"  and  "XI IF"  are  combined  into  "C1200FSK"  which  is  all  cases 

of  1200  baud  FSK  for  all  impairments.  This  class  has  been  assigned  the  plot 
symbol  "F".  In  a  similar  fashion,  Table  3-3  can  be  used  to  decode  any  of  the 
generic  name  information  of  Figure  5-27  and  the  remainder  of  the  OLPARS 
analysis. 

With  the  class  structure  defined,  the  next  step  is  to  perform  a  structure 
analysis  on  the  data.  Structure  analysis  is  performed  as  a  prelude  to 
classifier  design  in  order  to  verify  that  the  specified  tree  structure  (of  7 
olasses)  is  actually  present  in  the  data  as  defined  by  the  extracted  features. 
The  second  reason  for  structure  analysis  is  to  detect  the  presence  of  multi- 
modality  in  the  data.  Since  classifier  design  is  usually  dependent  on  the 
olass  statistics,  multi-modality  can  cause  gross  classifier  errors  depending 
on  the  particular  modality.  OLPARS  supports  a  variety  of  options  to  aid  the 
user  in  this  investigation,  A  subset  of  that  investigation  will  be  discussed 
here. 


After  the  means  and  covariance  matrices  are  computed  for  each  class,  the 
two  eigenvectors  corresponding  to  the  two  largest  eigenvalues  of  the 
oonvarienoe  matrix  are  used  as  a  projection  basis  into  a  two-dimensional 
space.  This  2-D  subspace  has  the  property  of  being  the  optimal  subspace  in  a 
least-squares  sense  (best  fit  to  the  current  data).  This  projection  (shown  in 
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Figure  5-28)  is  displayed  as  an  indication  of  the  class  clustering  and 
separation  between  different  classes  in  the  aoden  data  base. 


It  is  apparent  froa  the  projection  that  the  separation  between  the  seven 
classes  is  quite  good  based  on  the  features  used.  The  features  which  were 
used  in  the  eigenvector  projection  shown  in  Figure  5-28  are  lags  1,2, 3, 4,  and 
6  from  the  computed  autocorrelation  function  of  the  modem  data.  An  8K  sample 
window  size  was  applied  coupled  with  a  1-bit  quantization  of  the  data  in  order 
to  compute  the  autocorrelation  function. 

The  features  being  used  in  the  current  analysis  of  the  modem  data,  in 
OLPARS,  were  selected  from  a  total  of  30  possible  lags  of  the  autocorrelation 
functions  computed  on  the  MSA  data  base.  An  in-depth  analysis  was  performed 
on  all  of  the  features.  This  was  accomplished  by  using  the  discriminant 
measure  analysis  frame  in  OLPARS.  The  6  features  used  were  chosen  for  their 
ability  to  discriminate  between  the  different  modem  classes. 

Figure  5-29  illustrates  the  discriminant  analysis  output  from  OLPARS. 
The  analysis  is  broken  down  into  3  categories:  rank  by  logic  file,  rank  by 
class-alone  and  rank  by  class-pair.  The  ranking  of  the  features  is  based  on 
the  Fisher  direction  of  all  possible  class  pairs  in  N-space.  The  class-alone 
ranking  is  derived  by  summing  all  Fisher  direction  components  involving  the 
class  chosen.  The  rankings  of  all  features  are  listed  in  order  from  the 
largest  summed  Fisher  direction  magnitude  to  the  smallest.  The  class-pair 
ranking  is  computed  in  the  same  manner,  but  by  considering  only  the  Fisher 
directional  component  for  the  specified  class  pair.  Finally,  the  logic  file 
ranking  is  derived  from  the  sum  of  all  possible  class  pairs  listed  from  the 
largest  Fisher  direction  magnitude  to  the  smallest. 

From  Figure  5-29  it  is  seen  that  the  best  features  for  discriminating 
among  modems  are  1,4, 2, 6, and  3,  respectively.  The  best  features  for 
discriminating  class  C  or  the  CODEX  9800  modem  from  all  other  classes  are  1,4, 
2,3,  and  6,  while  features  4, 1,3*2,  and  5  are  the  best  choices  for 
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Figure  5-29:  Discriminative  Measure  Analysis  Features  (1,  2,  3,  4, 
5.  6,  7,  8)  1-Bit  Quantization  3K  Window  Size 
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discriminating  the  CODEX  LSI-9600  from  the  PARADYNE  LSI-9600,  based  on  the 
Fisher  directions.  While  Investigating  the  discriminant  analysis  for  best 
features,  other  classes  and  class-pairs  were  chosen  as  projection  planes  in 
the  analysis.  Features  1,2, 3. 4,  and  6  were  always  ranked  as  the  most 
discriminate  for  the  classes. 

Figure  5-30  illustrates  the  Fisher  pairwise  logic  evaluation  for  modem 
generic  type  separation  using  features  1,2, 3, 4,  and  6.  Shown  in  Figure  5-30 
is  a  confusion  matrix  for  classifying  each  class  correctly.  Choosing  class  A 
(HARRIS  5238  modem)  as  an  illustration,  it  is  shown  that  five  QAM  2667  baud 
feature  vectors  were  incorrectly  classified  as  QAM  2400  baud.  None  of  the  QAM 
2667  baud  data  was  rejected  thus  giving  a  correct  classification  of  91*23 
percent  and  an  incorrect  classification  of  8.77  percent.  An  overall 
performance  is  also  given  for  discriminating  between  all  classes  for  the  logic 
evaluation. 

The  Fisher  logic  evaluation  shows  that  98.8  percent  of  combined  design 
and  test  set  data  was  classified  correctly  and  1.20  percent  was  misdasslfled. 
An  extensive  analysis  was  performed  on  the  logic  evaluation  of  the  MSA  data 
base.  Figures  5-31  and  5-32  show  the  Eigenvector  projection  and  Fisher  logic 
result  for  the  same  features  (lags  1,2,3<4,  and  6)  but  using  a  3K  sample 
window  for  the  autocorrelation  estimation.  Note  the  increased  class  overlaps 
and  slight  decrease  in  overall  percentage  correct.  Analyzing  different 
combinations  of  window  sizes  and  quantization  of  the  data,  it  was  determined 
that  quantizing  the  data  to  1-bit  resolution  and  using  an  8K  window  size 
resulted  in  best  discrimination  and  modem  separation  for  features  1,2,3«4,  and 
6.  Using  a  1-bit  quantization  apparently  removes  the  amplitude  of  the  data  as 
a  feature  while  keeping  the  zero  crossing  information  intact.  The  resulting 
data  is  easier  to  separate  due  to  more  discrimination  in  the  computed  lags 
extracted  from  the  autocorrelation  function. 
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Figure  5-31:  Scatter  Plot  of  the  Eigenvector  Projection  for  the  5  Best  Lags  Computed  from 
the  3K  Point  Auto-Correlation  from  the  1-Bit  Quantized  Data 
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Figure  5-32:  Fisher  Logic  Evaluation  Based  on  Lags  1,  2,  3,  4,  6  of  the  1 -Bit  Quantized  Signal 

Auto-Correlation  Using  a  3K  Sample  Window 


A  similar  investigation  to  the  above  generic  typing  study  was  performed 
on  the  data  which  was  structured  by  baud  rate  (1200,  1600  ,  2400  ,  2667,  and 
4800)  and  then  structured  by  modulation  type  (QAM,  PSK,  PAM  and  FSK).  The 
respective  tree  displays  are  shown  in  Figures  5-33  and  5-34. 

The  results  of  the  three  investigations  are  summarized  in  Table  5-1.  The 
design  set  corresponds  to  ten  randomly  selected  vectors  of  each  lowest  sub¬ 
class  ,  for  example,  all  LENHJRT  26-C  modems  with  noise  added.  The  test  set 
Includes  all  of  the  Isolated  impairment  cases  together  with  the  oases  of 
combined  noise  and  phase  jitter  added. 

5.3  LATTICE  FILTERING 

Lattice  filtering  is  an  adaptive  linear  prediction  algorithm.  An 
adaptive  lattice  filter  attempts  to  set  prediction  coefficients  in  such  a  way 
that  the  difference  between  a  sampled  time  signal's  new  value  and  the  linear 
prediction  based  on  its  N  previous  values  is  minimized.  The  lattice  filter 
chooses  the  coefficients  in  a  time- symmetric  fashion,  so  that  it  predicts 
equally  well  the  sample  N+1  units  in  the  past,  given  the  most  recent  N 
samples. 

Another  interpretation  of  the  lattice  filter  is  that  it  models  a  signal's 
power  spectrum.  It  does  this  by  placing  zeros  in  the  z  plane  (the  sampled 
data  equivalent  of  the  frequency  domain)  so  that  if  the  zeros  were  replaced  by 
poles,  they  would  represent  a  good  model  of  the  signal's  power  spectrum.  The 
lattice  filter's  zeros  then  tend  to  cancel  the  peaks  in  the  signal's  power 
spectrum.  For  this  reason  the  adaptive  lattice  filter  is  a  spectrum  whltener. 

Figure  5-35  gives  a  schematic  diagram  of  a  lattice  filter.  The  signals 
along  the  upper  line  of  the  lattice  are  known  as  the  forward  residuals, 
because  each  represents  the  difference  between  the  actual  input  signal  and  a 
linear  prediction  of  it.  Similarly  the  signals  along  the  lower  line  of  the 
lattice  are  backward  residuals  with  each  representing  the  difference  between  a 
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figure  5-33:  Data  Tree  Illustration  for  the  Baud  Rate  Separability  Investigation 
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delayed  version  of  the  input  signal  8nd  an  estimate  of  it  based  only  on  the 
recent  past. 


In  the  figtre  each  box  marked  D  is  a  unit  delay.  Each  triangle 
represents  an  amplifier.  The  two  amplifiers  for  a  given  stage  have  the  same 
gain.  Circles  containing  pulses  are  summers.  Given  the  1th  forward  residual 
fi(t)  and  the  ith  backward  residual  b^  (t)  the  gain  ki+1  is  computed  according 
to 


ki-l 


-2Z  fi  (t)  b.  (t-1) 

Z  f?  (t)  +*b*  (t-1) 
t  i  t  i 


Then,  given  k^  +  1,  the  next  forward  and  backward  residuals  are  computed 
according  to 

fi+1(t)=fi(t)+ki+1bi(t_1) 

bi+1(t)ski+1fl(t)+bi(t~1)* 

In  its  capacity  as  a  spectrum  whltener,  the  adaptive  lattice  filter  is 
useful  for  the  present  effort  in  that  it  can  help  identify  the  type  of  modem, 
based  only  on  the  observation  of  a  noisy  signal  coming  from  it.  Whitening  the 
spectrum  of  a  signal  has  two  distinct  effects.  It  reduces  pulse  trains  into 
impulse  trains,  narrowing  and  sharpening  the  pulses.  It  also  decorrelates 
noise,  so  that  the  total  result  is  approximately  white  noise  with  short,  sharp 
pulses  mixed  in. 

The  itoitening  of  a  PS K  signal  then  tends  to  produce  a  series  of  pulses 
occurring  at  the  baud  rate.  Each  pulse  represents  the  transition  of  the 
carrier  signal  from  one  phase  to  another.  Some  pulses  will  be  positive, 
representing  one  polarity  of  phase  transition,  while  others  will  be  negative 
representing  the  opposite  polarity.  Then,  rectifying  the  lattice  forward 
residual  output  (for  PSK  input)  will  result  in  a  pulse  train  of  positive 


pulses  with  varying  amplitude,  with  the  repetition  rate  of  the  pulses  matching 
the  baud  rate. 

The  whitening  of  an  FSK  signal  tends  to  a  similar  result.  In  an  FSK 
signal,  a  tone  frequency  is  turned  on  and  off  while  a  tone  at  a  second 
frequency  is  turned  off  and  on.  This  is  done  so  that  one  or  the  other  tone, 
but  not  both,  is  on  at  all  times.  The  adaptive  lattice  filter  will  tend  to 
respond  (in  its  forward  residual)  with  impulses  at  the  points  where  the  tones 
are  turned  on  and  off  and  will  tend  to  block  out  the  tones  elsewhere.  Again, 
rectifying  will  result  in  a  series  of  pulses  at  the  baud  rate. 

Theoretical  study  indicates  that  two  stages  of  lattice  filter  suffice  for 
PSK,  and  four  stages  suffice  for  (two  tone)  FSK.  A  sinusoidal  signal  has  two 
poles  in  the  z  plane,  while  a  signal  composed  of  two  sinusoids  has  four  poles. 

A  lattice  filter  has  one  zero  per  stage,  so  it  can  fit  a  PSK  signal  with  two 

stages  and  an  FSK  signal  with  four. 

Experiments  were  performed  applying  lattice  filtering  to  the  modem 
signals  in  the  data  base.  Zero,  two  and  four  stages  of  lattice  filtering  were 
applied,  usually  followed  by  rectification  and  then  power  spectra  were  taken. 
Table  5-2  gives  the  results.  Line  height  refers  to  how  high  in  dB  a  line 
(spike)  in  the  power  spectrum  stood  out  above  its  neighbors.  All  but  two 
the  11  modems  tested  could  be  distinguished  by  generic  type.  Including  all  the 
lower  speed  modems.  While  the  initial  tests  were  run  on  low  noise  and 
distortion  signals,  later  tests  which  were  run  on  noisy  and  distorted  signals 
confirmed  the  results. 

Figure  5-36  through  5-43  illustrate  the  power  spectra  that  resulted  from 
the  processing.  For  example.  Figure  5-36  illustrates  the  result  of  four 
stages  of  lattice  filter  followed  by  rectification  and  then  taking  a  power 

spectrum ,  applied  to  the  CODEX  LSI-96  modem  signal.  The  spectrin  is  plotted 

on  a  logarithmic  scale  with  10dB  per  division.  A  spike  or  line  is  clearly 
visible  at  2400  Hz,  extending  about  5dB  above  its  surroundings.  In  some  of 
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ter  Spectrum  of  the  Rectified  Residual  from  Four-Stages 
Lattice  filtering  on  a  AOdB  Log  Scale 


ttimut/p 


Figure  5-38  Lenkurt  26-c,  Power  Spectrum  of  the  H*v  Signal  on  a  40  dB  Log  Scale 


Power  Spectrum  of  the  Rectified  Residual  from  Two-Stages 
of  Lattice  Filtering  on  a  40  dB  Log  Scale 


pectrum  of  the  Rectified  Residual  from  Four 
Lattice  Filtering  on  a  40  dB  Log  Scale 


the  plots  the  lines  in  the  spectrum  do  not  fall  exactly  at  their  nominal 
frequencies,  but  in  all  cases  they  are  within  16  Hz  of  their  nominal 
frequencies. 

Lattice  filtering  was  found  to  be  useful  in  identifying  modems  or  modem 
generic  type.  As  a  result,  it  is  useful  in  a  go/no-go  decision  process  to 
determine  if  a  received  signal  is  degraded  or  not  through  the  identification 
of  spectral  components  at  the  known  baud  and/or  carrier  associated  locations. 
Lattice  filtering  and  the  required  spectral  estimation  are  both 
computationally  complex  and  potentially  non-real izable  in  a  real-time  system. 

5.4  EYE  PATTERN  ANALYSIS 

The  hardware  eye-monitor  has  long  been  a  standard  industry  tool  in  the 
assessment  of  the  quality  of  analog  data  communications  links.  In  this 
section,  we  consider  the  aspects  of  a  software  implementation  of  an  eye- 
monitor  and  its  application  to  channel  quality  monitoring  of  digitized 
signals.  The  form  of  eye-monitor  investigated  was  the  I/Q  monitoring 
technique. 

The  received  VF  line  signal  is  assumed  to  be  of  the  form 

r(t)  =  Ret  n(t>  eJ  2ttc%\ 

where  f  is  the  carrier  frequency  and 

r(t)  =  r  (t)  -  j  r  (t) 

C  9 

is  the  complex  envelope  represented  in  terms  of  the  I/Q  components  r  (t)  and 
r  (t),  respectively.  The  I  and  Q  components  are  plotted  against  each  other  as 
shown  in  Figure  5-44  (a).  Bit  synchronization  with  the  input  data  stream  must 
be  achieved  and  used  to  sample  the  I/Q  components  at  the  mid-baud  Interval 
point.  This  is  the  point  where  a  decision  as  to  the  symbol  sent  would 


1 
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(a)  Oscilloscope  setup  for  eye 
pattern  monitoring:  the  I  channel 
provides  vertical  drive.  The  Q 
channel  horizontal  drive,  and  the 
system  data  clock  Z-axisj drive. 


(d)  X-Y  plot  showing  skewing  from 
phase  shifts. 


(b)  X-Y  plot  of  good  data  from  I  and  Q 
channels. 


(e)  X-Y  plot  showing  scattered 
points  due  to  transient  noise  and 
phase  hits. 


(c)  Plot  showing  effect  of  noisy  data. 

Figure  5-44:  Sample  I/Q  Plots  for  Common  Channel  Conditions 
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actually  be  made  In  the  receiver  portion  of  the  modem.  The  I/Q  plot  for  a 
good  data  channel  is  shown  in  5-44  (b).  The  effect  of  adding  noise  to  the 
channel  is  shown  in  Figure  5-44  (c).  Note  that  as  the  clusters  disperse,  the 
likelihood  of  a  sample  crossing  the  boundary  into  an  adjacent  region  and 
causing  a  symbol  decision  error  increases.  As  phase  jitter  is  Introduced  into 
the  system,  skewing  of  the  plot  occurs  as  illustrated  in  Figure  5-44  (d).  The 
effects  of  transient  noise  and  phase  hits  is  shown  in  Figure  5-44  (e).  In 
order  to  maintain  stability  in  the  I/Q  plot,  it  is  necessary  to  implement  some 
form  of  carrier  phase  tracking  to  prevent  random  phase  walk  or  mistuning  from 
rotating  the  I/Q  plot. 

For  the  purposes  of  this  investigation,  several  assumptions  were  made  to 
expedite  the  determination  of  the  effectiveness  of  this  approach  to 
performance  monitoring.  First,  in  lieu  of  carrier  phase  synchronization,  it 
was  assumed  that  the  carrier  was  constant  in  frequency  over  a  short  period  of 
approximately  60  ms.  With  this  assumption,  we  can  eliminate  the  requirement 
for  carrier  phase  tracking  if  we  know  the  precise  short  term  carrier 
frequency.  This  was  determined  using  the  Interactive  Digital  Receiver 
Simulator  (IDRS).  The  signal  was  demodulated  for  the  phase  waveform.  The 
complex  heterodyne  frequency  was  adjusted  to  show  a  phase  waveform  drift  of 
zero  over  the  time  window  being  monitored. 

The  problem  of  achieving  bit  synchronization  for  establishing  a  mid-baud 
interval  strobe  was  accomplished  in  the  following  manner.  The  signal  segment 
was  demodulated  for  the  FM  component  and  displayed  to  the  user.  Using  the 
cross-hair  input  of  the  graphic  terminal  (Tektronix  4014),  the  user  marked  the 
mid-baud  point  on  the  FM  waveform.  Since  the  baud  rate  was  known,  the 
software  could  sync  a  strobe  to  the  user-specified  point.  It  was  assumed  that 
the  baud  clock  of  the  modem  was  stable  over  the  60  ms  window  being  considered. 

Figure  5-45  shows  the  I/Q  plot  for  the  CODEX  LSI-48  Mode  B  modem  (DP5K, 
1200  baud).  As  the  signal  is  passed  through  the  4A  line  simulator,  the  I/Q 
plot  degrades  as  shown  in  Figure  5-46.  The  non-linear  phase  response  of  the 
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Figure  5-45  shows  the  1/Q  plot  for  the  CODEX  LSI-48  Mode  B  modem  (DPSK, 
1200  baud).  As  the  signal  is  passed  through  the  4A  line  simulator,  the  I/Q 
plot  degrades  as  shown  in  Figure  5-46.  The  non-linear  phase  response  of  the 
channel  characteristic  is  the  primary  cause  for  the  distortion.  The  other 
DPSK  1200  baud  signals  are  comparable. 

Figure  5-47  shows  the  I/Q  plot  for  the  CODEX  LSI-48  Mode  A  (DPSK,  1600 
baud)  for  the  no  4A  line  case.  Figure  5-48  shows  the  CODEX  LSI-48  Mode  C 
(QAM,  1600  baud)  also  for  the  no  4A  case.  For  both  of  these  modems,  the 
distortion  introduced  by  applying  the  4A  characteristic  was  sufficient  to 
distort  the  I/Q  plot  beyond  recognition.  Further,  none  of  the  higher  speed 
modem  signals  could  be  processed  to  present  a  stable  I/Q  plot  using  this 
technique. 

Several  conclusions  result.  The  requirement  for  some  form  of 
equalization  of  the  signal  as  the  baud  rate  increases  is  indicated.  The 
severity  of  the  group-delay  characteristic  is  sufficient  to  distort  the  phase 
plot.  Further,  the  stability  of  the  I/Q  plot  relies  on  carrier 
synchronization.  For  longer  term  monitoring,  carrier  phase  tracking  appears 
to  be  required.  Also  over  a  longer  sample  duration,  bit  synchronization  is 
required  as  a  strobe  indicating  when  to  sample  the  I/Q  components. 

The  next  section  of  this  report  proposes  an  alternative  method  of 
extracting  the  information  available  in  the  I/Q  plot. 

5.5  ORBITAL  VELOCITY 

One  possibility  for  signal  analysis  is  to  make  use  of  featires  extracted 
from  the  polar  coordinate  representations  of  the  in-phase/ quadrature  (I/Q) 
components  of  the  received  VF  line  signal.  For  want  of  a  better  name,  we  will 
refer  to  this  polar  plot  as  a  phase  trajectory.  In  the  present  subseotlon,  we 
describe  an  approach  to  performance  monitoring  or  evaluation  which  exploits  a 
feature  derived  from  the  phase  trajectory. 


Figure  5-45:  I/Q  Plot  of  the  CODEX  LSI-48  Mode  B  Modem  Signal  Without 

the  4A  Line  Simulator 
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The  received  V F  line  signal  is  assumed  to  be  of  the  fora 


r(t)  =  Re  (r  (t)e^  c  ) 


(1) 


where  ffi  is  the  carrier  frequency  and  ?f(t)sre(t)-Jrs(t)  is  the  complex 
envelope  represented  in  terms  of  the  I/Q  components  rc(t)  and  r9(t), 
respectively.  The  phase  trajectory  is  then  an  appropriately  normalized 
version  of  rs(t)  plotted  against  rc(t).  It  is  important  that  this  plot  be 
appropriately  normalized  to  avoid  coloration  by  an  unimportant  scale  factor. 
Specifically,  we  consider  plotting  the  quantity  (t)  sy(t)/o  ,  where 

a2  =  E{  |r  Ct)  |2}  =  E  {  r2  (t)}  +E(  rj  Ct)}  ,  C2) 

is  the  mean-square  value  of  the  complex  envelope  r(t).  is  also  the  mean- 
square  value  of  the  real  process  r(t).  The  quantity  9^  can  be  estimated  from 
the  received  line  signal  as  the  corresponding  sample  mean-square  value 
computed  over  some  finite  data  window. 

There  are  several  features  that  can  be  extracted  from  the  phase 
trajectory  that  should  prove  useful  in  performance  monitoring.  In  Figure  5-49 
we  illustrate  a  typical  phase  trajectory  associated  with  the  normalized 
complex  envelope  r'(t).  Here  the  point  r*(t)  on  the  trajectory  is 
parameterized  in  terms  of  its  polar  coordinates: 


a(t)  =  |  r  '  (t)  | 


C3) 


'V/ 


eCt)  =  arg  r  1  Ct) 


=-tan 


-Is 


(t) 


TtT 


C4) 


The  negative  sign  here  follows  from  the-  fact  that,  to  be  consistent  with 

r  1  |  I 

previous  work,  we  have  used  the  representation  r’(t)  *  rc(t)  -  JrsCt). 
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In  Figure  5-49  we  have  illustrated  only  one  circuit  of  the  complex  plane, 
whereas  the  phase  trajectory  associated  with  a  particular  observation  interval 
will  exhibit  many  such  circuits,  which  will  differ  from  each  other  for 
different  transmitted  data  sequences  and/or  channel  conditions.  The  phase 
trajectory  appears  to  be  useful  since  it  avoids  all  problems  of  phase 
synchronization  and/or  timing  recovery. 

Qne  of  the  features  readily  extracted  from  the  phase  trajectory  which  one 
might  consider  useful  in  performance  evaluation  is  the 
instantaneous  orbital  velocity  of  the  point  represented  by  r'(t)  as  it 
traverses  the  trajectory.  We  define  the  orbital  velocity  as 

s  (t)  =  ^2  A|7'(t)| 
dt 

1“  (t)+a~  (t)  O*'  (t)  ('  ) 

Here,  a(t)  is  merely  the  time  derivative  of  the  instantaneous  AM  waveform 
while 

6  it)  =  w(t) 


is  the  instantaneous  FM  waveform.  Both  of  these  quantities  are  easily 
extracted  from  the  normalized  complex  envelope  r'(t).  It  was  hypothesized  that 
second-order  statistics  extracted  from  s(t)  would  prove  useful  in  performance 
evaluation.  This  would  include  lagged  products  obtained  from  autocorrelation 
functions  and/or  spectral  properties  obtained  from  the  power  spectral  density 
function.  Note  that  if  the  phase  trajectory  is  along  an  arc  of  constant 
radius,  then  s(t)  i3  merely  proportional  to  the  magnitude  of  the  instantaneous 
radial  frequency.  The  use  of  orbital  velocity  in  this  sense  is  to  recognize 
that  motion  is  not  generally  along  an  arc  of  constant  radius  but  along  a 
highly  variable  and  sometimes  erratic  path.  This  might  be  due  to  filter 
transient  effects,  intersymbol  interference,  channel  noise,  etc. 


It  should  be  noted  that  In  making  use  of  orbital  velocity  we  have  assumed 
implicitly  that  any  residual  frequency  offsets  have  been  removed.  This 
requires  some  form  of  frequency  acquisition  and  tracking  scheme  which  should 
not  be  too  difficult  to  implement. 

The  parameter  s(t)  was  calculated  for  several  modems  of  the  database 
using  equation  (5).  Spectral  analysis  of  the  orbital  velocity  showed 
essentially  those  characteristics  seen  in  the  FM  power  spectrum. 

5.6  AM/FM  ANALYSIS 


We  have  seen  that  carrier  phase  instability  results  in  drift  in  the  modem 
signal's  phase  related  waveforms  and  I/Q  plots.  As  an  alternative,  the  FM  and 
AM  signal  components  were  thought  potentially  useful  features  in  the  analysis 
of  line  quality  due  to  the  immunity  to  mistiming  and  carrier  drift.  The 
mistuned  receiver  creates  a  drift  in  the  demodulated  phase  waveform.  The  FM 
component,  however,  being  the  derivative  of  the  phase  waveform,  shows  the 
mistuning  as  a  DC  offset  of  the  FM  signal  component.  The  AM  of  course  is 
determined  by  the  signal  envelope  and  is  not  effected  severely  by  mistuning. 

The  problem  of  carrier  drift  can  be  quite  severe.  The  specifications  for 
most  modems  require  carrier  accuracy  nominally  to  be  +3  Hz.  It  is  not 
uncommon  to  find  a  +5  Hz  range  in  practice. 

The  synchronous  demodulation  investigation  was  performed  using  the 

V*t 

Interactive  Digital  Receiver  Simulator  (IDRS).  The  result  of  the  demodulation 
of  the  WECO  207A2  modem  which  has  not  passed  through  the  4A  simulator  is  shown 
in  Figures  5-50  through  5-56.  A  frequency  shift  of  1800  Hz  was  used  with  a 
Chebyshev  1200  Hz  2-pole  low-pass  complex  filter.  No  smoothing  filter  was 
applied.  Figure  5-54,  showing  the  AM  power  spectrin  average  with  Hanning 
weighting,  clearly  shows  the  1200  Hz  baud  rate  component  and  harmonics.  The 
power  spectrum  average  of  the  rectified  FM  signal  shown  in  Figure  5-56  also 
displays  the  baud  rate. 
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e  5-51:  WECO  207A2  Demodulated  AH  Component  Haveform 


Figure  5-52:  WECO  207A2  Demodulated  Phase  Waveform 


5-56:  Power  Average  of  10  Lines  of  Rectified  WECO  207A2  FM 


It  ms  desired  to  investigate  the  effects  of  the  degrading  channel  on  the 
demodulated  AM  and  FM  components.  Using  the  spectral  average  of  the  AM  and 
the  rectified  FM  power  average  as  features,  the  baud  rate  components  were 
monitored  while  impairments  were  added.  In  the  case  of  the  DPSK  1200  baud 
modems,  no  significant  variation  in  the  AM  and  FM  spectrins  were  observed. 

With  the  higher  baud  rate  modems  it  becomes  difficult  to  perform 
demodulation  without  the  benefit  of  line  equalization.  No  unique  spectral 
components  were  observed  for  the  HARRIS  5238  or  PARADYNE  MP-96  modems.  The 
CODEX  LSI-96  AM  shows  a  2400  Hz  (baud  rate)  component  which  varies  as  a 
function  of  noise  added.  The  two  1600-baud  modems  (CODEX  LSI-48  Modes  A  and 

C)  both  showed  the  baud  rate  components  and  varied  only  slightly  with 
impairments  added. 

The  demodulation  analysis  of  the  FSK  modems,  however,  provided  some 
valuable  insight.  Although  no  stable  relationship  between  the  channel 
impairments  and  the  AM  and  FM  spectral  behavior  was  observed,  an  interesting 
case  was  discovered.  Figures  5-57  through  5-60  show  the  time  waveform, 
demodulated  phase  waveform,  FM  waveform  and  FM  power  spectra,  respectively, 
for  the  LENKURT  26-C  modem  with  phase  Jitter  added  at  10  Hz,  80°P/P.  Observe 
that  the  phase  waveform  (Figure  5-58)  shows  a  single  cycle  of  phase  shift  over 
a  period  of  1024  samples  (10  Hz).  The  amplitude  of  the  shift  is  estimated  at 
roughly  90°  P/P  from  the  vertical  axis.  The  FM  waveform  and  FM  spectrum  show 
no  such  contents,  however,  because  the  10  Hz  modulation  is  masked  by  the 
intentional  data  transitions  of  the  signal.  In  the  next  section  we  shall 
examine  a  technique  to  extract  this  unintentional  phase  jitter  modulation. 

In  summary,  the  spectral  analysis  of  the  AM  and  FM  derived  without 
carrier  phase  tracking  or  equalization  are  affected  by  channel  impairments, 
but  not  to  a  degree  which  is  readily  exploitable  or  exceptionally  stable. 


Figure  5-57:  LENKURT  26-C  Time  Waveform  Segment 


5.7  RESIDUAL  FM  ANALYSIS 


In  the  previous  section  we  observed  how  the  phase  jitter  impairment  is 
seen  in  the  phase  waveform.  Since  the  FM  is  the  derivative  of  the  phase,  the 
derivative  of  the  sinusoidal  variation  seen  in  the  phase  waveform  must  be 
present  in  the  FM  also  as  a  sinusoidal  component.  The  FM  recall  is  largely 
insensitive  to  carrier  drift  and  mistiming,  making  a  much  more  attractive 
basis  for  feature  extraction  than  the  phase  waveform. 

Figure  5-61  shows  the  histogram  of  the  demodulated  FM  component  of  the 
unimpaired  LENKURT  2 6-C  modem.  Since  most  time  is  spent  at  the  extremes  in 
the  FM  signal,  we  expect  the  "U"  shaped  profile.  The  means  of  the  histogram 
extremes,/^  and/^*  are  narked.  The  thresholds  T1  and  T2  are  defined  a  3/^/2 
and/^/g,  respectively. 

Based  on  the  histogram  information,  it  is  possible  to  effectively  remove 
the  intentional  data  transition.  Monitoring  the  FM  signal  amplitude,  we  can 
detect  whenever  threshold  Tj  is  crossed  and  subtract  out  the  expected  value  of 
the  data  transitions,/** 1 .  Similarly,  a  FM  transition  crossing  T2  is  detected 
by  thresholding  the  signal  at  T ^  and  subtracting  out  the  expected  value  of 
that  transition,//^  .  This  is  shown  in  Figure  5-62. 

Figure  5-63  shows  the  effect  of  removing  the  intentional  data  transitions 
on  the  LENKURT  26-C  sample  which  has  phase  Jitter  at  60  Hz,  80°P/P.  Low-pass 
filtering  this  residue  leaves  much  of  the  undistorted,  unintentional 
modulation  seen  as  five  periods/1024  sample  line  as  illustrated  in  Figure  5- 
64.  The  signal  is  resampled  by  4  and  a  high  resolution  power  spectrum 
generated  (Figure  5-65).  The  60  Hz  Jitter  component  is  enhanced  by  13dB  above 
the  noise  floor.  Figure  5-66  shows  the  spectrum  of  the  LENKURT  26-C  sample 
for  the  unimpaired  case  which  has  undergone  the  same  signal  processing  steps. 
Note  the  total  absence  of  60  Hz  spectral  components. 


FM  Signal 


Resultant  FM  Signal 

Figure  5-62  Illustration  Of  Data  Transition 
Removal  Via  Thresholding  By  T 
And  Subtracting  Out  The 
Anticipated  value, 


This  technique  was  applied  to  all  modem  types.  It  was  effective  in 
enhancing  the  Jitter  component  in  almost  all  cases.  The  results  for  the 
higher  speed  modems  (4800,  9600  and  1600)  usually  produced  a  measurable 
spectral  component  for  cases  of  50-80°  P/P  and  above.  Sample  FM  histograms 
used  for  the  determination  of  the  detection  thresholds  for  the  DPSK  1200-baud 
and  DPSK  1600-baud  are  shown  in  Figures  5-67  and  5-68.  For  the  DPSK  1200-baud 
modems,  the  detection  threshold  was  approximately  30°  P/P.  For  the  FSK 
modems,  phase  jitter  at  10°  P/P  could  generally  be  seen. 

A  shortcoming  of  the  algorithm,  however,  is  its  lack  of  sensitivity  in 
the  detection  of  the  10  Hz  phase  jitter  except  at  the  most  severe  cases  of 
60°-80°  P/P. 

5.8  ADAPTIVE  CHANNEL  MODELING 

In  adaptive  channel  modeling,  a  mathematical  model  having  a  number  of 
adjustable  parameters  is  formulated  to  represent  the  physical  transmission 
facility.  The  free  parameters  are  then  adjusted  by  adaptive  algorithms  to 
provide  a  close  fit  to  the  observed  signal.  Then  the  values  of  the  adapted 
free  parameters  describe  the  channel.  The  parameters  are  adapted  using 
decision  feedback.  That  is,  the  received  signal  is  demodulated,  the  data 
sequence  is  estimated,  and  the  result  is  fed  back  to  reduce  data-dependent 
effects  and  to  better  reveal  channel  effects. 

Appendices  A  and  B  present  studies  of  the  adaptive  modeling  concept. 
Appendix  A  derives  a  linear  filtering  model  of  the  channel,  while  Appendix  B 
derives  and  describes  a  nonlinear  model.  The  linear  filtering  model  is 
basically  an  adaptive  equalizer  combined  with  a  demodulator  and  remodulator. 
The  demodulator  determines  the  data  sequence,  and  the  remodulator  determines 
what  a  noise- free  undlstorted  version  of  the  signal  would  look  like.  Then  the 
adaptive  equalizer  estimates  what  linear  filtering  would  best  make  the 
original  signal  resemble  the  received  signal.  The  linear  filtering  is  modeled 
with  a  tapped  delay  line.  Once  the  linear  filtering  is  determined,  then  the 


difference  between  the  filtered  original  signal  and  the  received  signal  Is 
considered  additive  noise. 


The  adaptive  equalizer  can  be  replaced  with  a  non-linear  adaptive  device. 
A  method  for  doing  this  is  derived  In  Appendix  B.  The  method  Is  based  on  the 
Volterra  series  method  of  handling  nonlinearities.  The  resulting  adaptive 
nonlinear  device  is  shown  in  Figure  3  of  Appendix  B.  It  consists  of  two 
tapped  delay  lines.  One  is  used  to  model  linear  filtering  effects  of  the 
channel,  as  before.  The  other  has  its  output  squared  and  is  used  to  model 
quadratic  nonlinearities. 


The  adaptive  channel  models,  while  theoretically  interesting,  turn  out  to 
be  highly  complex  compared  to  the  requirements  of  the  present  application. 
Each  requires  a  model  of  both  the  modulator  and  demodulator  sections  of  each 
modem,  in  addition  to  adaptive  linear  or  nonlinear  devices.  One  possible 
method  of  implementation  might  then  be  to  use  a  bank  of  one  each  of  all  the 
modems  whose  signals  might  be  received.  An  alternative  would  be  to  model  all 
the  modems  using  software.  Due  to  complexity  constraints,  the  adaptive 
channel  models  were  not  implemented  in  the  present  effort. 

5.9  CHANNEL  SIGNAL  TO  NOISE  ESTIMATION 

The  ability  to  estimate  the  signal  to  noise  ratio  (SNR)  of  the  VF  line 
would  potentially  be  valuable  in  the  Phase  1  go/no-go  decision  as  well  as  in 
the  Phase  3  impairment  characterization.  Certainly  a  VF  line  which  exhibited 
an  exceptionally  poor  SNR  has  degraded  beyond  being  able  to  supply  acceptable 
service. 

Two  approaches  to  estimating  the  channel  SNR  were  considered:  spectral 
analysis  and  statistical  analysis.  The  spectral  approach  traditionally 
operates  on  a  sufficient  average  of  spectral  estimations  of  the  VF  line.  A 
region  of  the  spectral  average  is  identified  as  being  within  the  passband  of 
the  observed  process  but  not  as  being  used  by  the  system  for  data 
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transmission.  An  astlaata  of  this  idle  passband  region  power  is  foraed  by 
averaging  the  speotral  components  in  the  region.  A  similar  estimate  is 
performed  in  another  region  which  is  also  known  to  be  in  the  passband  of  the 
system  but  which  also  is  known  to  be  active  with  signal  transmission.  Since 
the  latter  estimation  contains  both  signal  and  noise,  the  estimate  of  the 
noise  power  is  subtracted  out  before  forming  the  SNR  as 

SNR  =  PS 

Pn 

where  Ps  and  Pn  represent  the  signal  power  and  noise  power  estimates, 
respectively. 


Examination  of  the  typical  signal  spectrum  from  Section  5.1  shows  that 
all  of  the  VF  line  bandwidth  is  utilized  by  the  signaling  strategy.  There  is 
typically  no  idle  band  in  the  VF  data  communications  system. 

This  suggests  the  requirement  for  statistical  estimation  of  the  channel 
SNR.  If  the  modem  input  is  assumed  to  be  a  random  sequence  (R(t))  then  over 
some  sufficiently  long  interval  At  we  expect  the  modem  transmitter  output 
(S(t) )  to  exhibit  the  same  output  amplitude  profile. 


As  noise  (N(t))  is  added,  the  composite  output  signal  C(t)  is  formed  as 
illustrated  in  Figure  5-69 .  We  now  define  the  peak  to  average  ratio  of  the 
composite  signal  over  the  N  sample  window  Tj  as 


pa  =  —2. 
PAj  /< 


c 


where is  defined  as  the  mean  of  the  full  wave  rectified  signal; 

Ac=  if  ?  ic(i)i 

i=l 


and  Cp  the  peak  value  of  C(t)  over  the  N  sample  window. 
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The  average  peak  to  average  ratio  over  all  N  sample  windows  of  the  signal 

sample  is  L 

^PA  =  “  £  PA. 

N  j-1  3 


over  the  L  available  N  sample  windows.  The  variance  of  the  rectified  peak  to 
average  ratio,  **PA,  is  given  as 


0PA  = 


N  E  A  -  PV: 


As  Gaussian  noise  is  added  to  the  system,  the  expected  value  of  the  rectified 

signal,  is  unchanged.  The  peak  value  of  C(t)  over  the  observation  window 
is  increased  by  the  peak  amplitude  of  N(t).  The  Increase  in  PAj  i9  therefore 
a  direct  function  of  the  amount  of  noise  added  to  the  channel.  The  average 
value  of  the  peak  to  average  ratio  (*pk  is  expected  to  be  a  stable  measure  of 
the  noise  present  on  the  channel.  This  is  Illustrated  in  Figure  5-70. 


If  the  N  sample  observation  window  is  not  of  infinite  duration,  then 
will  vary  in  proportion  to  the  probability  of  experiencing  a  maxlmun  value  of 
N(t)  over  the  observation  window.  In  the  absence  of  noise,  however,  l*pk  is 


quite  stable.  Therefore  we  expect  that  the  variance  of  H 
vary  in  proportion  to  the  amount  of  noise  present  on  the  line. 


,  will  also 


Let  us  now  consider  how  impairments  other  than  additive  Gaussian  noise 
will  affect  the  parameters'*^  and  °pA.  Referring  to  the  wireline  simulator 
functional  diagram  in  Figure  5-71 ,  it  is  observed  that  phase  jitter  can  be 
considered  to  be  low  frequency  FM  of  the  entire  signal  speotnjn.  The 
parameter  phase  jitter  frequency  controls  the  oscillator  frequency.  The  gain 
in  degrees  peak  to  peak  of  phase  jitter  controls  the  amount  of  FM  at  the 
specified  frequency.  FM  modulating  the  signal  directly  alters  the  signal  zero 
crossing  information.  The  signal  amplitude  is  affected  only  by  a  small  amount 
of  residual  AM  which  accompanies  the  intentional  FM  modulation. 


The  parameters  %  and  CFp^  are  dependent  only  upon  the  signal  amplitude 
Information.  Variations  In  the  zero  crossing  Information  of  the  signal  should 
not  affect  the  measurements.  Therefore,  the  signal  parameters  /ipA  and  apA 
are  largely  Immune  to  the  presence  of  phase  jitter  on  the  channel. 

The  effect  of  harmonic  distortion  on  AlpA  and  <7pA  Is  not  as  clear. 
Referring  again  to  the  functional  block  diagram  of  the  wireline  simulator 
shown  in  Figure  5-71,  we  observe  how  the  harmonic  distortion  Is  achieved.  The 
signal  S(t)  Is  passed  through  a  series  of  nonlinear  devices  which  raise  the 
input  signal  to  the  power  of  2, 3, J», 5,  etc.,  each  added  with  a  specified  gain 
to  the  input  signal.  The  resultant  composite  signal  C(t)  is  of  the  form 

4- 

C(t)  s  S(t)  +  aS2(t)  +  agS^U)  ♦  a^S^Ct)  ♦ 

where  aR  are  the  desired  gains  for  the  power  series  components. 

For  the  HSA  database  collection,  odd  harmonic  distortion  was  specified, 
the  effect  of  which  is  typical  of  the  application  of  a  soft-limiting  device. 
The  composite  output  signal  C(t)  is  of  the  form: 

C(t)  s  S(t)  ♦  a^U)  +  a4s5(t) 


Each  of  the  terms,  ajSN(t)  can  be  shown  to  contain  the  component  S(Nwt). 

The  odd  harmonic  distortion  which  occurs  during  soft  limiting  of  the 
signal  through  typical  telephone  exchanges  is  illustrated  in  Figure  5-7 2. 

Viewing  this  distortion  as  a  series  of  the  form 
C(t)  «  S(wQt)  ♦  a2S(3w0t)  +  a4S(5w0t) , 

the  resultant  waveform  is  being  squared  off  due  to  the  addition  of  the  higher 
in-phase  odd  harmonics.  The  more  severe  the  odd  harmonic  distortion,  the  more 
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constricted  the  resultant  waveform  and  the  lower  the  peak  to  average  ratio, 
^pA  as  well  as  the  variance 

From  the  above  discussion,  we  expect  that  the  effects  of  adding  noise, 
phase  Jitter  and  odd  harmonic  distortion  to  the  signal  will  affect  the 
parameters  /*pA  and  apA  as  follows: 

PARAMETERS  GAUSSIAN  NOISE  PHASE  JITTER  ODD  HARMONIC 

DISTORTION 

MpA  increase  no  change  decrease 

<rpA  increase  no  change  decrease 


An  extensive  experimental  Investigation  of  the  MSA  database  has  shown  this  to 
be  the  case. 

Software  was  developed  which  would  calculate,  for  each  waveform  file  of  a 
specified  list,  the  parameters  ^pA  and  CT  PA  for  a  specified  window  length.  A 
list  of  waveform  filenames  corresponding  to  all  samples  of  a  single  modem  was 
generated  for  each  of  the  11  modems.  Recalling  the  file  naming  convention 
described  in  Section  2,  the  list  of  all  LENKURT  26-C  modem  cuts,  for  example. 
Is  simply  all  files  beginning  with  the  character  "L".  The  modem  waveform  list 
included  all  impaired  cuts  for  that  modem  as  well  as  any  "clean"  cuts 
available.  The  list  of  modem  samp,  e  cuts  was  processed  to  produce  two  output 
files,  one  containing  the  estimated  values  of  MpA  for  each  cut,  and  the 
second  containing  the  <7pA  estimations.  A  separate  program  accepted  these 
files,  ordered  the  data  points  by  degree  of  impairment,  and  produced  an  output 
plot  for  each  parameter. 


Figure  5-73  and  Figure  5-74  show  these  plots  for  the  HUGHES  HC-276  2400 
bit/  second  DPSK  modem  using  a  window  length  of  256  samples.  The  left  axis 
corresponds  to  the  degree  of  impairment,  either  Gaussian  noise  or  harmonic 
distortion  added.  Harmonic  distortion  is  represented  on  the  plot  by  the 
symbol  'O'.  For  the  HC-276,  harmonic  distortion  was  added  over  a  range  from 
-30dB  to  a  worst  case  of  -4dB  in  8  steps.  Gaussian  noise  in  isolation  added 
to  the  signal  is  represented  by  the  symbol  "A",  and  is  present  from  -30dB  to 
-6dB  in  7  steps.  The  symbol  "^"  denotes  cuts  which  have  a  combination  of 
Gaussian  noise  and  phase  jitter  added.  Since  we  expect  that  the  phase  jitter 

will  not  affect  the  parameters  under  investigation,  the  cuts  containing  both 
impairments  are  plotted  only  by  the  Gaussian  noise  content.  Recall  that  any 
time  phase  jitter  occurs  in  combination  with  Gaussian  noise,  the  frequency  of 
the  phase  jitter  is  60  Hz. 

The  right  hand  axis  corresponds  to  the  isolated  phase  jitter  cases  of  all 
frequencies.  The  range  for  the  HC-276  is  from  10°P/P  to  80°P/P.  All 
frequencies  of  phase  jitter  are  denoted  by  the  same  symbol,  "0". 

Two  other  symbols  may  be  present.  First,  the  symbol  "•"  which  indicates  a 
cut  passed  through  the  wireline  simulator  4A  characteristic  but  with  no 
impairments  added.  Secondly,  the  symbol  "+"  which  is  a  direct  modem  sample 
by-passing  the  simulator. 

Examining  Figures  5-73  and  5-74  we  observe  the  "root"  of  the  function  on 
the  plot  base-line  corresponding  to  additive  Gaussian  noise  (AOI)  at  -30dB  and 
harmonic  distortion  (HD)  at  -30dB.  We  are  considering  a  channel  with  -30dB  of 
HD  or  AGN  to  be  completely  acceptable.  Note  that  as  HD  increases  on  the  line 
(the  symbol  □  ),  the  parameters  and  <7p^  decrease  in  a  regular  fashion. 
Similarly,  as  AGN  on  the  line  increases,  the  parameters  <y  pA  and^>A  increase 
regularly  regardless  of  the  frequency  or  amount  of  phase  jitter  (PJ)  combined 
with  the  noise. 
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Figure  5-74:  Peak  to  Average  Ratio  Variance  of  the  Rectified 
Signal  Waveform  for  the  Hughes  276  Modem 
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Figure  5-75  and  5-76  show  the  MpA  and  <7pA  plots  for  the  WECO  207*2  2*00 
baud,  DPSK  modem.  No  HD  was  collected  for  the  WECO  modem;  however,  the 
profile  of  the  plot  Is  as  anticipated.  The  increased  AGN  raises  the  0pA  and 
parameters.  The  PJ  cuts  are  within  a  small  range  corresponding  to  "good 
data" . 

The  CODEX  LSI-48  mode  B  2400  baud  DPSK  modem  is  shown  in  Figures  5-77  and 
5-78.  Again,  no  HD  was  collected.  The  PJ  and  AGN  points  exhibit  the  same 
characteristic  profile.  The  three  modems,  HC-276 ,  WECO  207A2,  and  CODEX  LSI- 

48  Node  B  shown  in  Figures  5-73  through  5-78  are  the  same  generic  modem  type; 
DPSK  2400  bits/  second,  1200  baud  with  a  1800  Hz  carrier.  We  would  hope  that 
the  characteristic  MpA  and  <xpA  functions  would  occupy  the  same  range  for 
each  modem.  Figure  5-79  and  5-80  show  the  three  modems  on  the  same  plot  for  a 
256  point  sample  window.  The  characteristic  profile  of  each  modem  is  slightly 
offset  from  the  rest.  Varying  the  length  of  the  window  used  in  the 
calculation  from  128  samples  to  1024  samples  did  not  improve  the  result.  It 
was  decided  that  some  subtle  variation  in  the  modem  signaling  format  or  the 
collection  process  caused  the  disagreement.  Software  was  developed  to 
normalize  the  input  signal  by  first  removing  any  bias  inadvertently  added  in 
the  collection  process,  and  second  by  normalizing  the  total  power  to  unity  in 
an  attempt  to  remove  any  discrepancies  in  collection  gain.  Some  improvement 
was  shown.  The  best  attempt  at  removing  the  offset  between  the  functions  is 
the  plot  of  the  normalized  signals  shown  in  Figures  5-79  and  5-80. 

Figures  5-81  and  5-82  show  the  MpA  and  0pA  plots  for  the  LENKURT  26-C 
modem  respectively.  Figures  5-83  and  5-84  show  the  /*pA  and  0pA  for  the 
MD-674  modem.  Both  modems  belong  to  the  same  generic  group:  FSK  modulation 
with  carriers  at  1200  Hz  and  2400  Hz  and  with  a  baud  rate  of  1200.  The  LENKURT 
26-C  plots  show  cases  of  harmonic  distortion  which  were  not  collected  for  the 
MD-674  modem. 
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Figure  5-77:  Peak  to  Average  Ratio  of  the  Rectified  Signal 
Waveform  for  the  CODEX  LSI -48  Mode  B  Modem 
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i  !  Figure  5-78:  Peak  to  Average  Ratio  Variance  of  the  Rectified 

|  Signal  Waveform  for  the  CODEX  LSI -48  Mode  B  Modem 
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Figure  5-79:  Peak  to  Average  Ratio  of  the  Rectified  Signal 
Waveform  for  the  HUGHES  HC-276,  WECO  207A2  and 
COOEX  LSI -48  Mode  B 
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Figure  5-81:  Peek  to  AveMge  Ratio  of  the  Rectified  Signal 
Waveform  for  the  LENKURT  26C  Modem 
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Figure  5-82:  Peak  to  Average  Ratio  Variance  of  the  Rectified 
Signal  Waveform  for  the  LENKURT  26C  Modem 
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Figure  5-83:  Peak  to  Average  Ratio  of  the  Rectified  Signal 
Waveform  for  the  MD-674  Modem 
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Figure  5-84:  Peak  to  Average  Ratio  Variance  of  the  Rectified 
Signal  Waveform  for  the  MD-674  Modem 
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The  cases  of  severe  harmonic  distortion  for  the  LENKURT  26-C  modem  are 
interesting.  Both  the  M pA  and  apA  functions  reverse  direction  for  cases  of 
severe  harmonic  distortion.  The  overall  characteristics  of  both  modems, 
however,  are  as  expected. 

Figures  5-85  and  5-86  show  the  two  FSK  modem  parameters,  shown  previously 
combined  on  the  same  plots.  Note  the  excessive  offset  between  the  "roots"  of 
the  functions.  There  are  several  possible  explanations  for  this.  First,  the 
duration  of  the  sample  MD-674  waveform  is  approximately  1.9  seconds,  as 
opposed  to  3.5  seconds  for  the  LENKURT  26-C  modem.  All  of  the  waveform  is 
used  in  the  calculation  of  the  statistics  so  the  longer  wave  will  produce  more 
stable  measurements.  Secondly,  the  HD-674  signal  has  been  re-sampled  down  to 
12.9  K  Hz  from  17.2  K  Hz.  Although  care  was  taken  in  the  resampling  process 
to  prevent  coloration  of  the  data,  the  possibility  exists.  Finally,  the 
histogram  of  the  time  domain  LENKURT  26-C  signal,  shown  in  Figure  5-87,  shows 
an  anomaly  in  the  symmetry  of  the  signal.  The  specific  cause  for  the  offsets 
of  the  roots  of  the  functions  of  Figures  5-85  and  5-86  has  not  been 
determined . 

In  actuality,  all  of  the  original  four  modems  suffered  a  bias  level  of 
approximately  20  units.  Further,  those  modems  had  a  dynamic  range 
corresponding  to  only  6  or  7  bits  out  of  the  12  available  from  the  A/D 
converter.  In  the  later  collection  of  7  modems  both  of  these  problems  were 
corrected . 

The  peak  to  average  related  statistics  for  the  CODEX  LSI-48  Node  A  modem 
(4800  bits/ sec,  DPSK,  1600  baud)  are  shown  in  Figures  5-88  and  5-89.  Although 
the  characteristic  plots  are  not  as  well  defined,  they  are  evident. 

The  same  statistics  for  the  CODEX  LSI-48  Mode  C  modem  (4800  bits/ sec, 
QAH,  1600  baud)  are  shown  in  Figures  5-90  and  5-91.  The  measurements  P  ^ 

a pA  are  unstable  for  this  modem. 
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Figure  5-86:  Peak  to  Average  Ratio  Variance  of  the  Rectified 

Signal  Waveform  for  the  LENKURT  26C  and  MO-674  Modems 
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Histogram  of  the  LENKURT  26-C  Modem  Time  Waveform 
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Figure  5-89:  Peak  to  Average  Ratio  Variance  of  the  Rectified 
Signal  Waveform  for  the  CODEX  LSI-48  Mode  A  Modem 
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Figure  5-90:  Peak  to  Average  Ratio  of  the  Rectified  Signal 
Waveform  for  the  C00EX  LSI-48  Mode  C  Modem 
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Figure  5-91:  Peak  to  Average  Ratio  Variance  of  the  Rectified 
Signal  Waveform  for  the  CODEX  LSI-48  Mode  C  Modem 
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Figures  5-92  and  5-93  show  the  plots  for  the  CODEX  LSI-96  nodes  (9600 
bits/ second  QAM,  2400  baud).  Note  that  even  at  9600  bits/ second  harnonic 
distortion  is  readily  detectable,  the  phase  jitter  region  is  bounded,  and  the 
average  of  the  noise  samples  tends  to  the  upper  right. 

The  PARADYNE  LSI-9600  modem  (9600  bits/second  PAM,  4800  baud)  is  shown  in 
Figures  5-94  and  5-95.  Harmonic  distortion  for  this  modem  is  detectable  at  18 
dB  and  greater.  Also,  the  phase  jitter  samples  are  well  behaved  while  the 
noise  samples  branch  off  in  the  l6dB  and  greater  region. 

The  plots  illustrating  the  PARADYNE  MP-96  (9600  bits/ second,  QAM,  2400 
baud)  and  HARRIS  5238  (16000  bits/ second,  QAM,  2667  baud)  modems  are  shown  in 
Figures  5-96  through  5-99.  For  these  two  modems,  as  with  the  CODEX  LSI-48 
Mode  C,  the  measurements  are  largely  ineffective. 

The  findings  of  this  investigation  appear  to  provide  valuable  insight 
into  the  estimation  of  channel  quality.  The  two  features  identified  are, 
first,  the  peak  to  average  ratio  of  the  rectified  signal  waveform,  and  second, 
the  variance  of  that  measurement  over  the  duration  of  the  signal  waveform 
sample.  The  effectiveness  of  these  measurements  in  a  demonstration  system 
shall  be  discussed  in  the  next  section. 
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Figure  5-92:  Peak  to  Average  Ratio  of  the  Rectified  Signal 
Waveform  for  the  COOEX  LSI-96  Modem 
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Figure  5-93:  Peak  to  Average  Ratio  Variance  of  the  Rectified 
Signal  Waveform  for  the  CODEX  LSI-96  Modem 
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Figure  5-94:  Peak  to  Average  Ratio  of  the  Rectified  Signal  Wave¬ 
form  for  the  PARADYNE  LSI-96  Modem 
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Figure  5-95: 


Peak  to  Average  Ratio  Variance 
Waveform  for  the  PARADYNE 


of  the  Rectified  Signal 
LSI -96  Modem 


5-125 


15.63  16.25  21.57  27.55  35.13  35.75  *4.35 

PHASE  JITTER  (oil  froqs,  P/P  dogroot) 


3.SI  9.SI  9.  SI  l.tt  9.54  1.S7 

SIGNAL  POWER  MEAN  VALUE  *11 

K  '  «  **  LI«  A  -  G0U33I0N  NO  I  SC  DOOCO 

«  -  CUOO  40  L1MC  y  -  MUSS  I M  NO  ISC  NITM  WMSf  JITTCO 

d  -  MMONNIC  OWOOTIM  000(0  o  *  190401(0  MM(  4JTU0.  DLL  rot 03 


Figure  5-97:  Peak  to  Average  Ratio  Variance  of  the  Rectified 
Signal  Waveform  for  the  PARADYNE  MP-96  Modem 
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Figure  5-98:  Peak  to  Average  Ratio  of  the  Rectified  Signal  Waveform 

for  the  HARRIS  5238  Modem 
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Figure  5-99:  Peak  to  Average  Ratio  Variance  of  the  Rectified 
Signal  Waveform  for  the  HARRIS  5238  Modem 
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6.  ALGORITHM  I  K»LE  STATION 

Based  on  the  investigations  described  previously,  several  algorithms  have 
been  found  which  appear  valuable  in  support  of  the  hierarchical  approach  to 
channel  quality  monitoring.  In  this  section  we  describe  the  implementation  of 
selected  algorithms  into  a  demonstration  system  capable  of  performing  gross 
channel  characterization,  generic  modem  traffic  classification,  and 
determination  of  the  type  and  degree  of  degradation  present  on  the  channel. 

The  system  is  implemented  in  a  modular  fashion,  a  single  task  performing 
a  specified  operation.  The  sequence  of  tasks  is  directed  by  the  Indirect 
Command  File  Processor  capability  under  the  RSX11M  operating  system. 

Figure  6-1  shows  the  block  diagram  of  Phase  1,  the  gross  channel 
characterization.  Phase  1  is  broken  up  into  two  distinct  operations.  Phase 
1A  is  the  peak/ average  ^PA)  and  peak/ average  variance  (°PA)  estimator  of  the 
rectified  input  signal.  For  all  the  modems  in  the  MSA  data  base,  these  two 
measures  fall  within  a  given  range.  This  range  information  is  stored  on  the 
disk  in  a  user- adjustable  file.  When  an  input  is  processed,  the  parameters 
1^4  and  °p£  are  estimated  for  a  256-sample  window  and  compared  against  the 
stored  range  Information.  If  both  measures  fall  within  the  given  range,  the 
channel  is  declared  "acceptable";  otherwise  it  is  rejected.  It  is  interesting 
to  note  that  test  tones,  voice,  and  most  noise  will  not  pass  because  of  the 
excessively  low  and  high  parameter  values. 

Phase  IB  is  not  used  in  the  go/ no  go  channel  decision;  rather,  it 
represents  a  collection  of  signal  processing  algorithms  which  have  been  shown 
to  be  generally  effective  in  enhancing  generic  parameter  information  (baud 
rate  and/or  carrier  frequency)  in  the  input  signal.  They  are  implemented  in 
an  optional  execution  mode  so  the  user  can  elect  to  receive  this  output  in 
graphic  form  for  his  analysis  purposes.  The  phase  IB  output  includes  the  raw 
and  rectified  signal  power  spectrun,  lattice  filter  2-stage  rectified  and  not 


6-1 


Imh* 


Figure  6-1 

PHASE  1  GO  /  NO-GO  DECISION  BLOCK  DIAGRAM 


rectified  power  spectrun  and  lattice  filter  4-stage  rectified  and  not 
rectified  power  spectrin.  These  six  spectrins  are  summarized  in  a  peaks 
analysis  report. 


Given  the  "acceptable"  channel  decision,  Phase  2  is  executed  to  classify 
the  generic  modem  type.  In  Phase  2  (Figure  6-2)  the  lags  1,  2,  3.  4,  and  6 
are  generated  for  the  autocorrelation  of  the  1-bit  quantized  input  signal.  An 
8000  or  3000  sample  correlation  window  is  used.  For  each  window  available 
from  the  input  signal ,  a  feature  vector  of  the  5  lags  is  generated  and  stored 
in  a  disk  file.  The  vector  file  is  then  processed  by  a  Fisher  discriminant 
pair  logic  classifier  against  3  logics,  the  outcome  of  which  are  the  BAUD, 
MODULATION,  and  combined  MODULATION/BAUD  class  indicators.  For  each  logic, 
the  majority  class  is  declared  the  result.  Given  that  either  of  the  first  two 
classifiers  agrees  with  the  combined  modulation/ baud  logic  result,  the  signal 
parameters  are  confirmed.  If  disagreement  in  the  logics  results,  the  decision 
is  that  the  signal  represents  either  a  bad  channel  or  an  unknown  modem.  In 
this  case,  no  further  analysis  can  be  done  on  this  signal.  In  an  operational 
system,  another  input  sample  could  be  requested  before  making  a  final 
determination. 

Given  a  decisive  output,  the  modulation  type,  baud  rate,  and  carrier 
frequency  (via  table  look-up)  are  available  for  further  signal  processing. 
The  classifiers  used  are  those  developed  previously  with  0LPARS  support.  The 
average  overall  percentage  correct  is  approximately  97+%.  Multiple  vectors 
and  multiple  classifiers  with  a  degree  of  certainty  weighting  in  the  review  of 
the  three  logic  results  provide  an  overall  generic  identification  of  high 
reliability. 

Phase  3  of  the  system  actually  performs  the  channel  quality  assessment 
(Figure  6-3 ).  Given  the  generic  ID  Information,  the  signal  is  FM  demodulated 
and  applied  to  a  data  transition  threshold  device  appropriate  for  that  modem 
type.  The  residue  is  200  Hz  low-pass  filtered  and  a  high  resolution  power 
spectrin  is  generated.  A  specialized  peaks  search  routine  deteots  phase 
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jitter  components  and  reports  them  to  the  user. 

In  addition,  the  signal  parameters  ^  and  crpA  are  used  in  a  look-up 
fashion  into  a  file  containing  a  template  of  those  parameters  for  this  modem. 
The  template  look-up  returns  the  type  of  impairment,  if  any,  (from  Gaussian 
noise  or  harmonic  distortion)  and  the  estimated  degree  of  the  impairment. 
Figure  6-4  shows  the  format  of  a  template  file,  either  for  UpA  or  0pA.  For 
any  calculated  measure  of  0r  °pA ,  there  is  only  one  corresponding  value 
in  the  template  file  which  directly  indicates  the  impairment  (HD,  GOOD  DATA, 
or  AGN )  and  the  amount  present.  The  results  of  the  demonstration  system 
implementation  are  discussed  in  the  next  section. 


6-6 


7.  RESULTS  OF  THE  DEMONSTRATION  SYSTEM 

The  demonstration  system  described  in  the  previous  section  was 
implemented  under  the  MSA  effort.  The  results  of  the  system  should  be 
considered  separately  for  each  of  the  three  phases  of  the  system  operation. 

Phase  1  of  the  system  performed  a  gross  channel -quality  analysis  on  the 
unknown  modem  signal.  The  peak-to-average  value  of  the  rectified  signal  and 
the  variance  of  that  measurement  over  approximately  two  seconds  of  input  data 
were  sufficient  to  make  the  decision.  The  range  of  these  parameters  was 
calculated  for  the  entire  database  and  was  used  as  an  acceptable  range  for  the 
channel.  Although  the  "classifier"  would  (by  design)  declare  all  of  the  MSA 
modems  as  "acceptable,"  it  does  reject  test  tones,  voice,  and  most  noise. 

Phase  2  of  the  system  performed  the  modem  generic  type  ID  by  applying 
three  Fisher  classifiers  to  sample  vectors  consisting  of  lags  1,2, 3, 4  and  6  of 
the  signal  auto-correlation  function  calculated  over  either  a  3000  or  8000 
sample  window  on  the  1-bit  quantized  input  signal.  The  individual  accuracy  of 
these  classifiers  was  described  in  Section  5.2.  The  three  classifiers, 
together  with  the  degree  of  certainty  logic  built  into  the  syst^<.  form  an 
exceptionally  accurate  ID  capability  for  the  signal  modulation  strategy  and 
baud  rate  Independent  of  the  type  and  degree  of  impairment  affecting  the 
signal . 

The  most  difficult  aspect  of  the  channel  monitoring  problem,  however,  is 
the  third  phase  of  the  hierarchy;  the  determination  of  the  type  and  degree  of 
impairment  degrading  the  line.  The  peak  to  average  ratio  of  the  rectified 
signal  and  the  variance  of  that  measurement  serve  as  a  basis  for  the 
determination  of  AGM  or  HD  on  the  line.  Although  this  technique  was  developed 
to  identify  impairment  at  the  generic  modem  level,  it  was  found  that  the 
variations  in  collection  parameters  and  the  post-collection  processing 
required  for  some  modems  introduced  sufficient  coloration  into  the  database  to 


mask  out  the  "fine-grain"  measurements  required  In  the  Impairment 
characterization.  It  was  shown  that  generically  Identical  modems  possess  the 
same  behavior  pattern  with  regarJ  to  the  measurements  extracted  but  that  they 
are  slightly  inconsistent  in  range.  For  this  reason,  a  fourth  classifier  was 
added  to  Phase  2  to  provide  the  actual  ID  of  the  modem.  The  ID  is  used  in 
selection  of  a  template  for  that  modem,  which  is  then  used  in  the 
determination  of  the  impairments  present. 

To  evaluate  the  effectiveness  of  the  AGN  and  HD  impairment  classifier, 
the  assumption  as  made  that  -25dB  and  less  AGN  present  or  -25dB  HD  present 
still  represented  a  "good"  quality  line.  The  decisions  of  the  two  classifiers 
were  treated  separately  when  counting  incorrect  decisions.  Table  7-1  shows 
the  ability  of  the  demonstration  system  to  detect  AGN  or  HD  on  the  in-service 
data  channel  to  the  accuracy  specified.  For  the  1200/2400  bit/ second  modems 
the  overall  ability  to  detect  AGN  or  HD  to  within  3dB  was  95.3%.  For  the  9600 
bit/second  CODEX  LSI-96  modem,  ACM  and  HD  were  identified  within  UdB  of  the 
actual  value  84.0%  of  the  time.  For  the  PARADYNE  LSI-96,  PARADYNE  NP-96  and 
HARRIS  5238  modems,  the  demonstration  system  was  not  able  to  accurately 
declare  the  impairments  present. 

The  system  was  successful  in  demonstrating  the  potential  usefulness  of 
the  phase  jitter  detection  algorithm.  The  detection  of  phase  Jitter  was 
accomplished  by  demodulating  for  the  FM  component,  threshold  detecting  for 
data  transitions  and  removing  them  by  subtracting  out  the  expected  transition 
value,  smoothing  the  residue  and  generating  a  high-resolution  power  spectrum. 
The  phase  jitter  component  is  increased  up  to  15dB.  Although  contract  time 
did  not  allow  the  mechanism  to  translate  the  spectral  result  to  actual  degrees 
of  jitter  developed,  a  high  resolution  peaks  analysis  routine  does  locate  the 
peaks  and  provide  all  necessary  amplitude  and  spectral  noise  measurements. 

Table  7-2  shows  the  detectable  levels  of  phase  jitter  for  the 
demonstration  system  even  in  the  presence  of  noise.  For  the  amount  specified, 
approximately  a  3dB  increase  over  the  noise  floor  for  the  tone  or  its 
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Table  7-1  (a)  Combined  AGN/HD  Classifiers  Overall  Percent  Correct 
by  Range  for  1200  and  2400  Bit/Second  Modems 


MODEM 

m 

+  2dB 

+  3dB 

Lenkurt  26-  C 

84.4% 

84.4% 

90.6% 

64 

MD-674 

87.0% 

92.6% 

98.2% 

54 

Hughes  HC-276 

96.6% 

96.6% 

98.3% 

58 

Code  LSI-48  Mode  B 

71.9% 

78.1% 

90.6% 

32 

Weco  207-A2 

85.4% 

91.7% 

97.9% 

48 

Overall 

86.3% 

89.5% 

95.3% 

256 

Table  7-1  (b)  Combined  AGN/HD  Classifiers  Overall  Percent  Correct 
by  Range  for  4800  Bit /Second  and  up  Modems 
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8.  CONCLUSIONS  AMD  RECO MjCNDATIONS 

Throughout  the  course  of  this  effort,  PAR  Technology  Corporation  has 
developed  and  demonstrated  algorithms  effective  in  the  monitoring  of  in- 
service  VF  data  lines.  A  hierarchical  approach  to  performance  monitoring  was 
presented  and  implemented  in  the  form  of  a  demonstration  software  system. 
Through  demonstration,  it  was  shown  that  reliable  performance  monitoring  can 
be  achieved  for  1200  and  2400  bits/second  modems.  Not  only  could  the  type  of 
impairment  be  recognized  fVom  among  Gaussian  noise,  odd  harmonic  distortion, 
and  phase  jitter,  but  the  degree  of  the  impairment  was  also  indicated.  Some 
success  was  also  shown  for  the  higher  speed  modems  although  the  accuracy 
decreased. 

The  results  of  this  effort  form  the  basis  for  further  Investigation  and 
show  the  practicality  of  an  eventual  hardware  implementation.  It  must  be 
realized,  however,  that  the  results  shown  under  this  effort  are  based  on  a 
limited  database.  First,  there  were  several  problems  with  the  database 
collection.  The  early  collection  effort  suffered  from  poor  dynamic  range  and 
mlsadjusted  A/D  converter.  The  overall  length  of  samples  was  not  consistent, 
nor  was  the  impairments  available  per  modem  sampled.  Finally,  two  different 
sample  rates  were  used:  12.8  KHz  and  17.2  KHz.  Although  the  latter  was  much 
more  desirable  for  the  "fine-grain"  Information  being  sought,  it  was  necessary 
to  down-sample  that  data  to  facilitate  using  the  entire  database  at  a  single 
sample  rate. 

The  next  step  toward  an  operational  performance  monitoring  capability  is 
the  evaluation  and  expansion  of  the  current  algorithms  on  a  more  complete 
database.  The  collection  must  be  performed  under  controlled  conditions, 
greatly  expanding  the  depth  currently  available  to  inplude  duplicates  when 
possible.  In  addition,  when  a  sample  modem  is  collected,  perhaps  the  signal 
could  be  put  over  an  actual  telephone  loop  back  circuit  and  then  collected. 
This  would  illustrate  immediately  the  effectiveness  of  algorithms  under  test 
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on  actual  data. 


In  review  of  the  algorithms  developed,  let  us  consider  each  phase 
briefly.  The  go/ no  go  channel  decision  performed  in  Phase  1  is  quite  accurate 
and  easily  computed.  The  generic  type  classifier  of  Phase  2  is  based  on  the 
1-bit  quantized  signal  autocorrelation  function.  This  is  simply  a  1-bit  A/D 
converter  in  a  hardware  implementation.  The  classifiers  are  computationally 
efficient  as  well  as  accurate.  The  Phase  3  AGN  and  HD  detection  as 
demonstrated  here  is  computationally  simple;  in  fact,  it  is  the  same 
computation  performed  in  Phase  1.  The  accuracy  can  probably  be  improved  upon 
simply  by  increasing  the  sample  rate  of  the  collection.  The  phase  Jitter 
detection  as  presented  here  is  rather  complex  and  computationally  difficult. 
A  more  suitable  approach  is  carrier  phase  tracking.  The  error  signal  from  the 
tracking  device  (vAiich  could  be  a  phase- locked-loop  or  PPL)  will  show  the 
phase  Jitter  component.  Although  more  difficult  to  implement,  the  impairment 
in  accuracy  should  well  be  worth  the  effort.  Further,  in  a  hardware 
implementation,  the  utilization  of  a  PPL  is  fairly  standard.  The  PPL  is 
equally  sensitive  to  low  frequency  Jitter  (lOhz)  as  well  as  higher  frequency 
Jitter .  It  is  not  difficult  to  see  how  this  technology  could  readily  be 
implemented  to  service  a  network  whose  population  was  4800  bit/ second  and 
slower  modems.  We  have  demonstrated  the  first  step.  The  expansion  and 
validation  of  the  current  algorithm  foundation  on  an  extensive  and  reliable 
database  is  the  bridge  toward  an  operational  capability. 
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Appendix  A 


Modeling  of  Voice  Frequency  digital 
Communication  Facilities  for  Performance  Monitoring  Purposes 


J.  V.  Modestino 

Electrical,  Computer  and  Systems  Engineering  Department 
Rensselaer  Polytechnic  Institute 
Troy,  Rev  York  12181 


I.  Introduction: 

There  are  several  techniques  available  for  performance  monitoring  of  voice 
frequency  (v?)  digital  cosauni cations  channels  as  discussed  in  [l].  One  such 
technique  is  the  use  of  adaptive  channel  modeling  techniques  as  described,  for 
example ,  in  [2j-[3]  for  use  on  the  mlcrovave  line-cf-sight  (LOS)  channel. 
Specifically,  a  finite-parameter  model  of  the  physical  transmission  facility 
is  constructed.  The  free  parameters  of  this  model  are  then  adjusted  by  an  appro¬ 
priate  adaptive  algorithm  to  provide  a  dose  fit  to  a  given  channel  chosen 

for  performance  evaluation.  In  the  present  note  ve  rill  attempt  to  develop 
an  appropriate  finite-parameter  model  for  the  digital  VF  channel  where  the 
use  of  the  word  channel  in  this  context  is  intended  to  include  both  the  modem 
transmitter  as  well  as  actual  transmission  facilities.  That  is,  we  are  inter¬ 
ested  in  degradations  introduced  by  the  modem  itself  as  well  as  the  physical 
propagation  channel.  The  issue  of  how  the  free  parameters  are  to  be  adapted 
to  a  particular  channel  are  described  only  briefly  here. 

H.  Preliminaries.: 

In  what  follows  ve  make  extensive  use  of  complex  narrowband  representations 
of  all  signal  and  noise  processes.  Specifically,  the  received  signal  to  be  used 
for  adaptive  channel  estimation  is  given  by 

r(t)  *  Re{r(t)eJ2irfct}  ,  (1) 


A- 1 


(2) 


where  Ms  u  isssed  to  ova  carrier  frecuency  la  S z.  and 

C  " 

?(t)  ■  rc(t)-Jrf(t)  , 

is  the  ccsalex  envelope  expressed  in  terns  of  lowpass  waveforms  re(t)  and 
r((t)  representing  the  inohase  and  Quadrature  (X/Q)  components ,  respectively. 

The  eocplex  enrelepe  f(t)  is  furthermore  as  coned  of  the  fora 

«(t)  ■  *0(t)  +  fi(t)  ,  (3) 

where  SQ(t)  is  the  complex  envelope  of  the  received  steal  component  and  fi(t) 
is  a  complex  noise  process  representing  additive  channel  noise  as  veil  as 
possible  adjacent  channel  interference  effects.  The  received  signal  compon¬ 
ent  lQ(t )  is  assured  modeled  as  the  complex  response  of  the  channel  to  the 
transmitted  sir-*--  !(t)  at  its  input.  For  purposes  of  the  present  discussion, 
the  channel  is  modeled  as  a  linear  time-invariant  narrcvband  filter  with  complex 
impulse  response  function  he(t)  as  indicated  in  Fig.  1.  In  reality,  the  channel 
might  veil  he  nonlinear  and/or  time  varying,  nevertheless,  the  channel  can  et 
least  be  expected  to  he  time-invariant  over  intervals  that  are  long  relative 
to  typical  measurement  intervals.  The  question  of  channel  nonlinearities  is 
samevhat  more  difficult  to  address.  Eovever,  since  most  of  the  modulation 
strategies  of  interest  possess  constant  envelopes  the  presence  of  narrowband 
nonlinearities  is  not  expected  to  have  a  large  influence  on  the  performance  of 
the  digital  data  link.  At  any  rate  the  question  of  channel  nonlinearities 
remains  to  he  mere  fully  explored  in  ongoing  work. 

Ve  assume  that  the  transmitted  signal  !(t)  is  normalized  so  that 

Eg  -  |T,|,(t)|zdt,  (U) 

0 

represents  the  signal  energy  where  T#  is  the  baud  interval.  Furthermore,  v* 
restrict  attention  to  digital  signaling  formats  such  that  l(t)  has  the 


Sere,  hB(t)  Is  s  complex  waveform  vanishing  outside  the  interval  [o,T#]  and 
{S^}  is  a  complex  sequence  which  depends  upon  the  data  to  he  communicated  in 
a  manner  specific  to  the  modulation  strategy  employed.  For  example,  the  case 
Sj*  ±1  and  ha(t)*UQ(t)  vith  Ug(t)  the  real  baseband  pulse  waveform 


u0(t> 


;  0<t<T 


elsevhers 


corresponds  to  conventional  coherent  BPSK  aodulation.  Similarly,  arbitrary 
KPSK  can  he  handled  hy  appropriate  choice  of  {S^}.  The  quantities  t  and  6  in 
(5)  are  random  timing  epoch  and  phase  offset,  respectively.  Generally  these 
will  he  functions  of  time.  For  the  time  being  we  will  assume  that  both  t  and 
6  are  known  at  the  receiver  to  arbitrary  accuracy.  Without  loss  of  generality 
then  we  asaume  both  x  and  6  are  zero  in  (5). 

The  received  signal  component  can  then  he  represented  over  the  measurement 


interval  t-UT^HTj]  as 

l«(t) 


ft.  H 


l0(t)7  tJ  iI_IT  5i  (T) 

where  h(t)^hB*he(t)  is  the  convolution  of  the  complex  modulator  vaveform  h^Ct) 

end.  the  complex  impulse  response  he(t)  of  the  channel.  The  overall  impulse 

•  • 

response  h(t)  is,  of  course,  not  known  since  the  channel  response  h  (t)  is  in 

c 

general  unknown.  Nevertheless,  It  may  he  possible  to  adaptively  estimate  this 
quantity  which  should  provide  useful  information  on  the  present  state  of  the 
channel.  A  scheme  for  accomplishing  this  is  illustrated  in  Fig.  2.  Eere  the 
complex  envelope  f(t)  is  passed  through  a  filter  with  a  complex  impulse  response 


h„(t)  and  th tz  stalled  it  the  output  la  synchronise  -1th  the  baud  rite.  The 

tiding  inf omit ion  vill  presucably  cose  from  i  hit  synchroniser  vhich  rest 

he  replieited  is  pirt  of  the  modem.  The  stapled  output  is  then  compared 

* 

to  m  estlaated  sequence  {f^}  obtained  by  applying  the  estiaited  data  sequence 

A 

{ei>  to  a  finite  pinaster  model  of  the  chtanel.  The  estimated  lata  sequence 
* 

{51>  vill  be  assusied  to  cone  froa  the  aodea  output  vhich  is  replieited  is  part 
of  the  performance  monitoring  unit.  Zt  is  shove  is  Fig.  2  as  being  obtained  it 
the  output  of  an  appropriately  defined  decision  device.  Ve  assuae  that  the 
error  rate  is  reasonably  small  so  that  the  occasional  errors  do  not  have  a  pro¬ 
found  effect  on  overall  system  performance.  Assuring  that  an  appropriate  finite. 

* 

parameter  model  of  the  channel  is  available,  the  expected  output  sequence  {f^} 
is  coopered  to  the  actual  sequence  {f^}  to  generate  the  error  sequence 
An  appropriate  adaptive  algorithm  is  then  employed  to  estimate  and  trad  changes 
in  the  channel  paraaetars.  Clearly,  the  reason  for  including  the  receiving  filter 
in  Fig.  2  is  to  obtain  a  SHR  advantage  in  the  comparison  of  expected  and  actual 
channel  outputs. 

In  order  to  determine  an  appropriate  finite-parameter  model  for  the  channel, 
note  that  the  sampled  output  sequence  is  given  hy 

f,  -  I  f<v)h  (t-t)dT  ;  i-0,il,±2,...,  (8) 

J.  r  t«iTf 

or  in  normalised  form 


fi  V  V?T  Jl*  **  hi-k +  fti  ’  (9) 

vbere  is  a  complex  zero-aean  Gaussian  sequence  vhose  I/Q  components  are 


each  of  unit  variance. 


<>*•**» 


‘  ' 


3  i 


(10) 


is  the  equivalent  rectangular  bandwidth  of  the  receiving  filter  is  Hz, 

finally 


h^*  j"h(t-iTs)hr(J?i-t)4t,  i,J-0,±lt...  (11) 


The  sailed  output  sequence  can  be  revritten  in  the  fora 


2Z  [2 JT~  H 

yT  ^Vi-k^i 

fc*i 


(12) 


vhere  ve  assise  ':,«1  vhich  is  merely  a  conservation  of  energy  constraint.  The 
second  tern  in  (12)  is  a  distortion  term  resulting  from  ISI  produced  by  co- 
channel  distortions  on  the  channel. 

Tran  the  preceding  it  is  dear  that  a  reasonable  model  for  the  channel  is 
the  tapped  delay  line  (TEL)  illustrated  in  Tig.  3-  Specifically,  ve  assume  the 
channel  can  be  modeled  as  a  TEL  vith  a  total  of  2L+1  taps  and  corresponding 
tap  veights  S4 ,  i*0,±l,i2,...,±L.  In  the  next  section  ve  describe  hov  the  tap 

veights  nay  be  adjusted  to  estimate  and  trach  the  instantaneous  state  of  the 


channel. 


HI.  Adaptive  Adjustment  Procedure: 

Define  the  performance  criterion  as 

J(*)  »  Klfj-fJ2}  (13) 

m 

vhere  £  *(g_j),8_Ij+1,. .  .,J0, . . .  ,4^)  represents  the  tap  gain  vector  of  dimension 
(2L+1)  and  5^  represents  the  TDL  filter  output  according  to 


L 


<  it) 


t  Sere  8  (j«0is  the _ system  transfer  function  of  the  receiving  filter  vith 
lapulsS  response  h„(t). 


where  for  siapiicixy  we  assume  tb*.t  ,  i.e.,  so  decision  errors  Bade. 

The  performance  criterion,  which  is  to  he  nirinized  through  choice  of  £, 
can  he  expressed  u 


•  M|»t|»>-aw  l_L^nh-i 
♦I  I  l,*I ■ 

J«-L  k«-L  3  “  3  1  * 


» 


(15) 


This  last  result  esc  he  rewritten  in  rector  matrix  fora  according  to 

j(£)-s(  |?i 1 2)-i+g"VCi-4~1^+S^-2"1i]  »  (16) 

where  £  is  a  (2L+l)x(2L+l)  matrix  with  (h,i)  component 

\,»*  ««L»W  *  n,k  •  <»> 

and  *  is  a  (2L+1)  rector  with  k’th  component 

*k  -  ;  k-0,±l . ±L  (18) 

The  "t"  in  (16)  indicates  complex  conjugate  transpose.  It  is  easily  seen  that 
the  optimum  tap  connection  rector  is  that  which  minimizes  the  quadratic  fora 
in  (lS)  with  the  result 

ion  *  i'1*  •  <»> 

This  solution  can  he  determined  hy  an  iteratire  gradient  procedure  provided  the 
gradient  of  the  functional  J(^)  can  he  determined.  Specifically,  we  propose  use 
of  the  steetest  descent  procedure  where  £q  is  an  initial  guess  and  the  successive 
estimates  are  determined  according  to 

“  k— 1 »  i*0»lf-» 

with  h  the  step  size  at  each  iteration  and  the  gradient  2Si)  -*  determined 
according  to 


(20) 


!<*)  * 


£ 

3£ 


t 


(21) 


cr  equivalently, 

I(&)»2[§t-*1  »'-22{*i£1)  .  (22) 

*  m 

where  is  the  error  signal  and  8^* (6^,8^^. . .  ,8^,. .» ,8^) 

is  the  transnitted  data  sequence. 

The  difficulties  with  implementing  the  steepest  descent  algorithm  as 
described  by  (20)  is  that  the  gradient  in  (22)  requires  an  ensemble  average. 

As  an  alternate  ve  propose  use  of  the  "noisy*'  gradient  estimate 

1^)  =  -21^8*]  ,  (23) 

vhich  is  the  typical  approach  used  in  practice. 

IV.  Buggery  end  Conclusions; 

We  have  described  one  approach  to  modeling  and  adaptive  estimation  sf  the 
characteristics  of  VF  communication  lines.  Several  things  are  required  to  com¬ 
plete  this  vork.  First,  techniques  should  be  developed  for  estimating  perform¬ 
ance  analytically  in  terns  of  the  estlaated  parameter  vector  &.  Also  pattern 
recognition  techniques  should  be  explored  for  classifying  the  degradation  using 
the  vector  £  e*  a  pattern  rector.  Sene  experimental  results  should  be  obtained 
using  this  approach.  Finally,  ve  should  investigate  aore  accurate  modeling 
approaches  to  include  nonlinear  effects  in  a  finite-parameter  channel  model. 
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I.  Introduction : 

In  previous  work  Cl] ,  an  adaptive  approach  to  performance  monitoring 
of  digital  communication  systems  operating  on  voice  frequency  (7?)  lines 
vas  described.  This  approach,  vhich  vas  an  initial  attempt  at  formulating 
and  solving  the  problem,  is  limited  in  a  number  of  vays.  Most  important 
among  these  limitations  are  the  fact  that  the  modulation  format  must  be 
knovn  a  priori  and  a  replica  of  the  actual  receiver  modem  be  constructed, 
and  the  fact  that  only  linear  channel  distortions  can  be  treated  in  this 
say.  It  is  of  some  interest  to  devise  alternative  techniques,  adaptive  or 
otherwise,  which  remove  either  or  both  of  these  limitations.  In  the  pre¬ 
sent  note  we  describe  one  such  technique  which  at  least  allows  modeling  of 
nonlinear  distortion  effects. 

H.  Preliminaries : 

The  Volterra  series  expansion  can  be  used  to  represent  the  output  of  a 
nonlinear  system  with  memory  in  powers  of  the  input.  This  fact  has  been  very 
useful  in  the  analysis  of  the  input /output  behavior  of  nonlinear  baseband 
systems  as  dtscribad,  for  example,  in  The  extension  to  narrowband 

signals  described  in  terms  of  inphase  and  quadrature  (I/Q)  components  is 
relatively  straightforward.  In  particular,  suppose  that  x(t)  and  y(t)  are 
narrowband  signals  represantlag  the  input  and  output  respectively  of  a 


narrowband,  al though  possibly  nonlinear,  system  as  illustrated  la  Fig.  1. 
Hers  vs  assume  that  x(t)  and  y(t)  can  bs  represented  in  tbs  form 

x(t)  -  ^ReUCt^V)  ,  (la) 

and 

y(t)  *  /S  Se{?(t)eJ“ct>  , 

vitb  tii-  the  carrier  frequency  in  radians /second  and  2(t),  y(t)  tbs  corrsa- 
c 

ponding  complex  envelopes.  The  input  /output  behavior  of  this  system  is 
assumed  described  by  the  complex  Volterra  series 


?(t) 


•r 


dt  b  (t-Tl,t-T2, 

an 


,t-r TTPrtO. 

a  i*l  1 


(2) 


vbere  the  b  (t. ,t_ . t  )  are  called  Volterra  kernels.  For  a  linear  system 

n  i  2  n  — 

all  terms  except  the  first  vanish  and  ve  have 

tit)  »  £  h^(t— t )X(t)dt  ,  (3) 

•  «• 

in  vhicb  case  h^(t)-h(t)  is  simply  the  complex  envelope  of  the  impulse  re¬ 
sponse  h(t)  of  the  linear  narrowband  filter.  The  impulse  response  h(t)  then 
has  the  representation 

h(t)  •  2  Re{h(t)e^wct>  .  (U) 

Another  special  case  occurs  when  the  system  is  nonlinear  but  memoryless. 

It  follows  that  the  Volterra  kernels  are  described  by 
“  ,a, 

h  (t-r i,t-T2,...ft-T  )  *  all  «(t*r  )  (5) 

n  a  1-1  i 

where  i(*)  is  the  delta  function.  From  (2)  then  we  have 

tit)  -  l  af^t)  ,  (6) 

n-1 


which  is  merely  a  representation  of  the  output  as  a  power  series  expansion 


of  the  input  complex  envelope.  One  cia  conclave  other  instantaneous  or 


2aro-*aaor7  noalinearicies  vhich  distort  the  amplitude  sad  phase  ia  decid¬ 
edly  diff treat  vays.  These  import sat  special  eases  vill  he  treated  ia  later 
wort. 

Finally,  consider  the  a-fold  Fourier  transform 


Hn(«l dt2***|" dt^Ct^.tg,. . .  ,ta)emp{-j(a»i,t1+W2t2 


+...  (a  t))  .  (7) 

U  IX 


Clearly  H^(<u)  is  the  familiar  complex  transfer  function  of  a  narrovband  linear 
system  given  by 

H^Cee)  -  j  h^tje^dt  .  (8) 


By  analogy,  the  quantity  H  .  ,«a)  vill  be  called  the  a'th-order  com¬ 

plex  Tolterra  transfer  function. 

HT.  Application  to  Adaptive  p«i-»nT-«»nce  Monitoring: 

As  in  [l]  we  shall  assume  the  received  signal  to  be  of  the  form 

r<t)  «  G  Re{f(t)eJ“ct>  ,  (9) 

with 


f(t)  »  iQ(t)+a(t)  ,  do) 

share  lQ(t)  is  the  complex  envelope  of  the  received  signal  component  end 
fi(t)  Is  a  complex  noise  process  representing  additive  channel  noise  as  veil 
as  possible  adjacent  channel  interference  effects.  The  received  signal 
component  ig(t)  is  assumed  modeled  as  the  complex  response  of  the  channel 
to  the  transmit*—*  a(t)  at  its  input.  We  assume  the  transmitted  signal 

is  of  the  fora 

,(t)  9J\l  fiiVt-iVT),J0  »  (U) 
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7«* 


with  hm(t)  4  complex  waveform  vanishing  outside  the  interval  [o.T^],  and 
(c. }  a  complex  sequence  which  depends  upon  the  data  to  be  communicated  in  a 
Banner  specific  to  the  sodulatlon  strategy  employed.  The  quantities  t  end 
9  in  (11)  are  random  timing  epoch  and  phase  offset,  respectively.  As  in 
(l],  we  will  assume  for  the  time  being  that  both  r  and  9  are  known  to  arbit- 
rary  accuracy  at  the  receiver.  Without  loss  of  generality  then  we  assume 
that  both  t  and  9  in  (ll)  are  zero. 

The  channel  behavior  will  be  modeled  in  terms  of  a  truncated  complex 
Volterra  series.  Specifically,  we  assume  that 


SQ(t) 


)TT  s(o,  (12) 

i«l  1 


for  same  choice  of  h  (t),  n*l,2,...,K.  For  ease  of  exposition  we  take  the 

c,n 

case  X«2.  in  what  follows .  The  case  of  finite  K '  follows  in  similar  fashion 
except  the  algebra  is  a  bit  tedious.  Thus 


IQ(t)  *  J~  hc(t-t)s(T)dT+|*  |*  hc(t-T,t-v)l(t)S(\,)drdv  , 


(13) 


where  we  have  suppressed  the  dependence  upon^  n  in  writing  h  ( • )  and  h  ( • ,  • ) 

c  c 

for  h  ,(•)  and  h  „(•,•),  respectively.  The  overall  TO*  channel  model  is 
c»x  c,z 

then  as  illustrated  in  Fig.  2. 

The  received  signal  component  can  then  be  represented  over  the  measure- 
meat  interval  [-ST^.CJ+ljT^]  as 


H  E  IS  V 

(t)  .fii  l  a  h(t-iT  )  *=2-  l  f  5  S  h(t-iT  ,t-JT  )  , 
1  *  s  i--a  *  u  i— h  i— u  1  a  *  * 


(1M 


t  This  should  cause  no  confusion  since  the  linear  and  quadratic  filtering 
functions  are  readily  identified  by  vhether  it  is  a  function  of  one  or 
two  arguments. 


•  A  •  • 

where  h(t)  *  ha*b,  ( t )  is  the  convolution  of  the  complex  modulator  vavefom 


b^(t)  and  the  complex  Impulse  rttpout  h^(t)  of  the  channel.  Similarly ,  the 

quantity  h(t,s)  is  the  two-dimensional  (2-D)  convolution  h(t,s)^  h  **h  (t,s) 

s  c 


defined  by 


h(t,s)  *  |*  |*  bB(t-t)hc(t,v)ha(s-v)4Tdv  . 


(15) 


It  should  be  obvious  hcv  additions!  high-erder  convolution  products  arise  if 
the  Voltaire  sarlas  of  (12)  is  truncatad  for  sons  value  3>2.  In  the  ease 

9  m 

under  consideration  neither  h(  • )  or  h(  • ,  • )  is  known  and  aust  be  adaptively 
estimated  from  the  output  of  the  7?  channel  - . 

As  in  [l],  we  assume  that  the.  channel  output  f(t)  is  passed  through  a 
narrowband  filter  with  complex  impulse  response  h^(t)  and  sampled  in  synchron¬ 


ism  with  the  baud  rate  to  generate  the  sequence 

;  i-O.tl^,..., 

t«lT 


■  f*  f(t)h^(t-T)dT 


(16) 


or  in  normalized  form 

2E_  H 


/  2Es  H  .  /  22  J  B 

h  -75^7  J.,  Vi-k  Vs^7  *!_,  tL,  Vi-k-i-t't  **i 


(17) 


Bare,  {n^}  is  a  complex  zero-mean  Gaussian  sequence  whose  I/Q  components  are 
each  of  unit  variance, 

B  ^  2t  |*  ,  (18) 

is  the  equivalent  rectangular  bandwidth  of  the  receiving  filter  in  Hz, 

fct.j  *  I"  h(t-JTa)hr(iTs-t)dt,  . . . ,  (19) 

and  finally, 
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*i-i,i-Al  J  *<*-«i»*-*VV1Vt)4fc  i  (20) 


It  should  be  noted  at  this  point  that  ve  have  written  these  last  two  express¬ 
ions  only  in  texas  of  i-J  In  the  ease  of  (19)  and  the  double  index  i-3c,  i-t 
in  the  ease  of  (20).  Justification  for  this,  which  is  not  immediately  appar¬ 
ent,  is  relatively  straightforward. 

Trcaa  the  preceding,  it  is  clear  that  a  reasonable  finite-parameter  model 


of  the  channel  is  of  the  foxa 

L 


fl '  jit,  *  *L  JL  5i-A,  A-i 


(a) 


The  first  term  has  an  interpretation  as  a  conventional  tapped  delay  line  (TBL) 
as  in  Cl] ,  while  the  second  texa  cn  the  right-hand  side  of  (a)  represents  a 
quadratic  operation  on  the  data  sequence  which  is  a  bit  more  difficult  to 
Interpret. 

One  interpretation  to  be  given  to  the  quadratic  operation  in  (a)  is  to 


express  it  in  texas  of  vector  matrix  notation  according  to 
u  L 


Xl  tL  <J1’  *  St?  ’ 


(22) 


vh€?6  2^  *  (2^^ •  *  •  i2^* •  •  *  9  S  ^  *  (ZL+l)x(ZL+l) 

•• 

matrix  with  (k,l)  component  A,  k,i«0,xl,*2, . . . ,  L.  We  will  assume  that  g 
is  positive  definite  and  symmetric  in  which  case  it  can  be  factored -in  the 
foxa  g  ■  Jpg  where  g  is  an  upper  triangular  matrix.  It  follows  that 

L  L  •**.*■«  ■  ^  • 


(23) 


where  *  g  5^  .  In  fact,  under  the  preceding  assumptions  we  have 


(2k) 
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jjjjich  is  la  general  _  *  tine-varying  filtering  operation.  For  example,  with 
we  hare 

L  . 

!i  '  i  Vt'i-i  (a5> 

*  l  Vi-t 

i»0  4  1  » 

viiar •  E.  ^  b  ..  We  will  assume  that  the  filtering  operation  represented  by 

*  Ot^ 

the  matrix  g  oea  he  approximated  by  the  input /output  relationship  in  (25) . 

*»  • 

Thi«  Till  he  the  ease,  of  course,  if  h^  A*hQ  which  we  assume  to  he 
approximately  true.  The  sequence  (t£>  then  has  an  interpretation  as  the 
impulse  response  sequence  associated  with  a  linear,  causal  and  time-* invariant 
finite- impulse  response  (FIR)  digital  filter  ,  which,  can  again  he  realised  as 
a  TDL-  Indaed,  the  overall  channel  model  can  he  realized  as  a  parallel  com¬ 
bination  of  two  TDI/s  with  a  squaring  device  and  an  accumulator  associatad 
with  the  quadratic  branch,  A  possible  realization  of  this  structure  is 
illustrated  in  Fig.  3.  As  indicated,  there  are  a  total  of  (3L+2)  coefficients 
to  he  specified*  (2L+i)  coefficients  and  (L+l)  coefficients  (E^}^>Q. 

In  the  next  section  we  propose  an  adaptive  procedure  for  estimating  these 
coefficients  on  the  basis  of  past  transmissions. 

17.  Adaptive  Adjustment  Procedure: 


As  Indicated  In  the  preceding  section,  the  output  of  the  nonlinear  channal 
modal  can  he  expressed  in  the  fora"*" 


t  L  ft  T2 

h  “  l  Vi-j  *  I  l  Vi-j-k  m  i  . 

1  J-_L  J  1  J  j— L  k*0  *  1  J  *  1  1 


t  We  are  implicitly  assuming^  decision  directed  schema  d.«d  as  discussed 
previously  in  Cl]  •  Also  (f«'>  represents  the  output  of  theHiaear  branch 
while  (f  I ' }  represents  the  Sutput  of  the  quadratic  branch  as  illustrated 


Sow  define  the  performance  measure 


(2T) 


m  j 

where  |  ■  . if)  *ad  i  *  (iQ,  ^,...,1^)  represent  the 


tap  gala  rectors  in  Fig.  3  of  dimension  (2L+1)  tad  ( L+l ) ,  respectively.  This 
performance  measure  is  to  be  minimized  by  simultaneous  choice  of  g  sad  i,  and 


esa  be  rewritten  as 


J(g,§)  •  SClffJ 2)-2He{£  f^fje  » 

J«-L 


+2He 


IL'do  Jo  *« 

0L 

'll  l 

■  I  I  *-*  I  S(?*5  ,5,  1  a>  ’ 

[m*l  a-0  *  a  j—t  1  i’jHa  i  J~Q 


-2He 

L  L 


(28) 


1 1  I  VS  I  I  VS  Z  Z  S(ii.jHl!l.H!U-a!U.T}  • 

m-0  a«0  “u-C  r-0  J*-L  k--L  1  J  1  J  1  M  1  * 
While  this  expression  looks  quite  complicated,  it  is  shown  in  Appendix  A  that 


it  caa  be  simplified  somewhat  and  rewritten  in  rector  matrix  notation  accord¬ 
ing  to  ^ 

J(  g ,  I  )-E(  |  t±  i  A)  ]  +i~L( i-X(  a)  ] 

♦ci-i’1(  j-i(a) )  )  i 

-23S(tr(TtiT]  HtrCll^U)  >  (29) 

Zaeh  of  the  terms  in  this  last  expression  are  defined  explicitly  in  Appendix  B, 
to  which  the  reader  is  referred  to  for  details.  Furthermore,  it  follows  easily 

t  Sere  the  dagger  "t"  indicates  complex  conjugate  transpose. 


that  the  optimum  choice  for  |  and  i  a re  found  as  aolutiona  to  the  simultaneous 


equations 

*  W(|.*)  - 

•  o  ,  (30*) 


sad 


♦ 


w(f.*) 

» 

URe{£(i)*}  -  0 


(30b) 


la  general,  these  simultaneous  nonlinear  equations  are  difficult  to  solve, 
ffeverthaless ,  the  solution  can  be  determined  by  an  iterative  gradient  pro¬ 
cedure.  Specifically,  we  propose  use  of  a  steepest  descent  procedure  where 
|q  aad  ig  are  initial  guesses  and  the  successive  estimate*  are  deterained 
according  to 


It+i  *  51-h!2J(Sf!i)  *»  wAf*  » 

(31a) 

5i+i*ii"il2*J’(ti’-i)  ;  * 

(31b) 

where  h  is  the  step  sise  at  each  iteration  aad  the  padienta  are  detesalned 
from  (30). 

To  provide  some  physical  interpretation  to  this  procedure,  consider  the 
gradient  of  (30a),  i.e. , 

y(|,I)  »  2Se{g|-(J-X (a)))-  -auKKiJ*})  (32) 

•  T  / 

where  is  the  error  signal  and,  as  previously,  *'ai>L,5i+L-l’ 

. ..,ai,...,ai_t(+i,ai>>E()  is  the  trannitted  data  sequence.  Similarly,  we  have 

TjJ’Cf  .i)  -  J*R2{2  '•tSi»  > 
where  now  b.  is  an  (L+l) -vector  with  Jt'th  element 


(33) 


L 

^  .5i-jai-j-ic  *’  ,  (3M 

J*-L 

with  the  sequence  of  outputs  of  the  upper  tapped  delay  line  la  Fig.  3 

with  components  generated  according  to  (25) ,  i.e., 

L 

aj  ,  -  I  lna.  ;  J-0,±l,h2,  . . . ,±L  .  (35) 

i-J  n-0  a  1-J-a 

The  difficulties  with  implementing  the  steepest  descent  algoritlsa  as 
described  above  is  that  the  gradients  in  (32)  and  (33)  require  ensemble 
averages.  As  an  alternate  we  propose  'use  of  the  "noisy”  gradient  estimates 


VgJ(8,4)  =  -2He{«ia*>  , 

(36a) 

y(|,I)  2  -i*Re{«J*}  , 

(36b) 

which  is  the  approach  usually  adopted  in  practice. 

V.  Sunnary  and  Conclusions; 

Ve  have  described  an  approach  to  modeling  and  adaptive  estimation  of  the 
characteristics  of  VF  communication  lines,  subject  to  nonlinear  distortion 
effects.  As  noted  previously  in  [  1  ],  it  remains  to  develop  analytical  tech¬ 
niques  for  estimating  error  rate  performance  In  terns  of  the  estimated  para¬ 
meter-  vectors  £  and  i  .  Also  to  be  investigated  are  pattern  recognition 
techniques  for  classifying  the  degradation  using  the  vectors  £  and  A  as 
parameter  vectors.  Finally,  we  note  that  the  approach  described  here  is  limit¬ 
ed  to  linear  modulation  strategies .  It  remains  to  determine  how-  adaptive 
estimation,  techniques  can  be  used  with  nonlinear  modulation  formats . 


I 
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Appendix  A 

Vector-Matrix  7oraula.ti.on  of  Perforaance  Criterion 


Connider  tie  perforaance  criterion  J(f,i)  u  given  by  (28)  of  tie  text. 


This  can  be  written  as 


J(|,a)  »  +  J2(-)*  (. 

*****  i*Peada  upon  boti  g  and  i  and  is  obtained  as  tie  sun  of  tie 

first  four  teras  in  (28),  i.e., 

*  /  < 

J.(g,i)»E{ !f  i2>-2Se'  J  g. S{r*c.  .)  •  • 

_  J»-l  J  1 


ft  L  L  L 

+2Sa  ’  I  84  l  l  I  s(a.  .«*  .a*  .  } 

j«-L  -a-0  n-0  a  “k--!,  i”J 


(A-l) 


L  L 


(A-2) 


SlatHarly,  tie  term  J2(i)  depends  only  upon  I  and  consists  of  tie  last  two 
teras  of  (28),  i.e.. 


L  L 


J2(I)»  -  2Re-  l  l  II  [  E{f*d  5 
2  [n-0  n-0  a  aji-L  1  i'jHa  i_J- 


*  Jo  n^O^Jo  ioWjLt  J.LS{ai-JHa5i-j-u?i-k-n5*-k-T}-(A-3: 

We  concentrate  upon  each  of  tiese  teras  separately. 

After  some  algebra  it  is  easily  shown  that 

^(f.IJ-EC (fj 2>-(j-i(§)]tj"1[i-X(i) ]*{|-S’1(5-X(C) )]tg(g-^*1(i-X(l) ) ]  .  (A-10 
^ » •  •  •  *^q»*  •  •  »^)  is  a  ( 21+1 )  —vector  with  k'th  ccaponeat 


B-ll 


Similarly,  $  is  a  (2L+l)x(2I-*-l)  matrix  with  (k,l)  component 


(A -6) 


Finally,  X^CaW.X^X^+i . .  ,Xg,. . .  fX^)  is  a  (2L+1)- vector  with  Jc'th  com¬ 
ponent 

K  *  1,  ivJ  <A-T> 

»"0  n*0  1«-L 

and,  we  have  written  X(£)  as  an  explicit  function  of  £  to  emphasize  this  dep 
ence  as  indicated  by  (A-7) .  It  should  he  noted  in  particular  that  the  gradie* 


with  respect  to  g  is  given  by 


Similarly, 


sj  (|,£)  .  _  ,  -  . 

71Jl(|,£)  S  -  2Se{g[£-£  1(<j-X(£))]> 


'  -V- 


(A-8) 


*  URe{$(f)£>  , 

with  |(|)  an  (L+l)x(L+l)  matrix  depending  upon  |  with  (nx,n)  sleaent 


L  L 


■•a1*’  *  '  *U(*>  ' 


(A-9) 


(A-9) 


(A-10)| 


Sow  from  the  definition  of  Jg(£)  111  (A**3)  we  can  rewrite  this  in  the 


r2(l)“  -  25e{tr(yaiT]>  *  trU*+A(i))  . 

Sere  f  is  a  (L+l)x(I.+i)  matrix  with  (m,n)  element 

£  r  Y  - 

V» "  jl_t  *  *»,*  i  . «■ 


(A— 11) 


(A-12) 


t  Here  it  is  easier  to  go  bach  to  the  definition  of  given  in  (A-2). 


Finally,  the  matrix  £(i)  is  an  (I+l)x(L+l)  matrix  depending  explicitly  upon 
i  vitk  (a, a)  element 


x* 

n,a 


(A-13) 


It  follovs  after  some  elementary  algebra  that 

3J-(I) 

ltfzU)  -  -jy-  -  -  !*Re{»J}+  URe{A(i)a>  . 


U-iU) 
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Appendix  C 


Appendix  C 


Modem  Signature  Analysis  Demonstration  System  Users  Guide 


C.1  Overview 

The  MSA  demonstration  system  is  implemented  as  modular  tasks  linked 
together  with  the  RSX11M  Indirect  Command  File  processor.  This  facilitates 
the  removal  of  certain  modules  for  other  applications  without  code  modifica¬ 
tion. 

On  the  DICEF  PDP  11/40  system,  the  MSA  directory  is  [3,3J.  Command  files 
exist  together  with  the  sources  in  that  directory  so  that  a  complete  rebuild 
of  the  software  system  can  be  effected  with  the  command; 

«  MSABlD 

The  rebuild  would  only  be  necessary  in  cases  of  operating  system  modifi¬ 
cation. 

In  order  to  operate  the  MSA  demonstration  system,  the  user  must  be  in  the 
MSA  directory  and  have  the  desired  input  waveform  files  in  the  same  directory. 
The  waveform  is  expected  to  be  in  a  fixed  (IDRS)  format  file  which  has  the 
following  attributes: 

-  unformatted 

-  direct  access  organization 

-  512  bytes  per  record 


Reoord  1  of  the  file  is  a  header  containing  the  information  shown  in  Table  C- 
1.  This  data  format  is  harmonious  with  most  A/D  oolleotion  systems.  Data  is 


C-1 


__ . 


stored  from  record  2  on,  1  sample  per  two-byte  Integer  word.  Also  It  Is 
assumed  that  the  waveform  Is  sampled  at  12.8  Kc. 

C.2  Operation  of  the  MSA  Demonstration  System 

With  the  12.8Hz  sampled  waveform  in  the  MSA  directory,  [3.3J.  the  MSA 
command  file  is  invoked  with  the  command 

9  MODEM. 

The  interaction  is  now  explained. 

•DO  YOU  WISH  TO  SKIP  PHASE  IB?  [Y/N]: 

Optionally  the  generation  of  the  raw  and  rectified  signal  spectra  and  the 
power  spectra  for  the  2-  and  4-stage  lattice  filtered  signals,  both  rectified 
and  not,  can  be  omitted  from  the  experiment.  Since  none  of  this  output  is 
used  in  the  phase  1  go/no  go  decision  these  calculations  are  not  required. 

•AUTOCORRELATION  BUFFER  LENGTH  (3K  OR  8K): 

Here  the  user  specifies  the  correlation  length  as  "3K"  or  N8K”  points. 
Note  that  more  vectors  result  for  the  ”3K”  decision. 

•AUTOMATED  [Y/N]: 

The  system  can  process  in  batch  mode  if  the  data  is  located  on  magtape 
"mounted"  on  the  tape  drive  and  a  list  of  the  desired  magtape  files  is  con¬ 
tained  in  "WAVENAME. LST"  in  the  format  illustrated  by  the  command  file  prompt. 


•WAVEFORM  FILE  TO  BE  PROCESSED? 


This  requests  the  name  of  the  input  signal  to  process,  assuming  that 
automated  mode  was  not  selected.  This  is  the  end  of  the  user  interaction. 

Output  from  a  sample  run  is  shown  in  Figures  C-1  through  C-1 1  for  the 
CODEX  LSI-48  Mode  B  modem  with  60  Hz  phase  jitter  at  64°  P/P  added.  The  out¬ 
put  is  self-documenting  with  the  possible  exception  of  the  peaks  detector  sum¬ 
maries  which  are  described  in  the  next  section. 

C.3  MSA  Demonstration  System  Program  Reference  Manual 

The  Command  File  Structure 

The  'MSA  Command  File'  is  divided  into  3  main  phases: 

1.  Gross  Channel  Quality  Analysis 

2.  Modem  Generic  Type  ID 

3.  Channel  Impairment  Characterization 

Each  of  these  phases  is  further  subdivided  into  smaller  tasks,  all  linked 
by  the  RSX11M  Indirect  Command  File  Processor.  The  structure  is  detailed  in 
the  following  section. 

C.3. 1  Phase  1:  Go/No  Go  Decision 

C.3.1'1  Phase  1A:  Gross  Channel  frailty  Estimation 

The  Go/No  Go  decision  of  the  channel  is  performed  using  time  domain 
statistics: 

-  Rectified  Peak/Average  Signal  Ratio 


Rectified  Signal  Peak/ Aver age  Ratio  Variance 


Programs: 

NORMALIZE  -  normalize  inputs  the  specified  raw  IDRS  formatted  waveform, 
converts  to  R*4  data  as  required,  removes  any  residual  bias 
caused  by  the  A/D  Converter,  and  normalizes  the  waveform  to 
unity  power.  The  output  waveform  has  the  extension  ".NRM". 

EVLSTS  -  accepts  the  (".NRM")  normalized  waveform  and  computes  the  rec¬ 

tified  signal  peak  to  average  ratio  and  rectified  signal  peak 
to  average  ratio  variance. 

A  reference  file  "EVLSTS. DA T"  contains  the  user-settable 
acceptable  limits  for  each  of  these  (2)  tests.  The  output  from 
the  program  is  the  PHASE  LA  summary  and  a  flag  passed  back  to 
the  Indirect  Command  File  Processor  signaling  either  'CHANNEL 
ACCEPTED'  or  ’CHANNEL  REJECTED'. 

C.3. 1.2  Phase  IB:  Signal  Spectral  Enhancement 

Phase  IB  is  a  collection  of  signal  processing  algorithms  deemed  promising 
in  enhancing  certain  spectral  components,  such  as  baud  or  carrier,  to  a  point 
where  they  are  suitable  for  classification  features  for  generic  typing  or 
impairment  characterization.  This  phase  consists  of  (3)  major  algorithms 
linked  as  shown  in  Figure  C-1. 

Programs: 

RECT  -  accepts  an  input  waveform  and  creates  a  full  wave  rectified 

output  waveform.  Input  and  output  filenames  are  specified. 


C-A 


MSLFLT 


is  the  multi-stage  lattice  filter.  Input  and  output  waveforms 
are  specified  with  the  desired  number  of  filter  stages. 


PVRSPT  -  is  the  spectral  estimator.  The  input  waveform  is  specified  and 

successive  1024-point  real-to-oomplex  FFT's  are  performed  and 
averaged  to  a  maximum  of  10.  The  output  is  in  waveform  format. 
512  points  (bins)  long. 

PEAKDT  -  is  the  peaks  detector  which  operates  on  the  spectrum  created  by 

PWRSPT.  Only  specified  frequency  ranges  of  the  spectrum  are 
considered  as  noted  on  the  output  listing.  Each  of  the  6  spec¬ 
trins  generated  in  Phase  IB  are  processed  concurrently,  and  the 
results  are  consolidated  onto  8  single  output  summary  page. 

A  sample  cell  from  the  peaks  detector  has  the  following  format 
and  fields: 

M  2170.  Hz 
C  0.  Hz 
65.db, .  Hz 

where: 

"H"  s  maximum  frequency  in  band 
"C"  s  center  of  valid  spectral  spike 
"§tt  s  minimum  clear  pulse  height  above  noise  in  order 
that  a  spike  be  declared  valid 
" s  BW  of  the  valid  component  "$"dB  down. 

NOTE:  If  no  acceptable  spike  wes  found, 

"C"  and  values  will  be  zero. 
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C.3.2  Phase  2:  Modem  Generic  Type  Classifier 


Five  lags  of  the  autocorrelation  (1,2, 3, 4  and  6)  of  the  1-bit  quantized 
signal  are  generated  for  the  user-specified  window  of  3K  or  8K  samples.  All 
of  the  normalized  signal  is  processed  and  the  output  vector  file  ( ,VEC)  con¬ 
tains  one  vector  of  5  lag  measurements  each  for  each  non-overlapping  window 
available  in  the  input  waveform. 

The  vectors  are  run  through  four  Fisher  classifiers: 

-  Modulation  Type 

-  Baud  Rate 

-  Combined  Modulation  Type  -  Baud  Rate 

-  Modem  ID 


ACORR  -  generates  specified  lags  of  the  autocorrelation  of  the  normalized 
input  waveform.  Lag  0  is  used  for  normalizing  the  vector  (other 
lags).  The  program  creates  a  ".VEC"  file  suitable  for  input 
directly  into  the  Fisher  classifier  (or  OLPARS).  The  length  of 
each  correlation  window  is  user-specified  at  3000  or  8000  samples. 


FHREXE  -  is  the  Fisher  Discriminate  pair  classifier.  Input  files  are 
".VEC"  vector  files  and  a  summary  is  produced  together  with  flags 
signaling  the  .AT  processor  as  to  the  outcome.  The  logic  library 
is  user  specified. 


PHASE2  -  Phase2  is  the  executive  program  which  compiles  summary  and  effects 
final  Phase  2  logic  decision  of  modulation  and  baud. 

If  the  baud  rate  or  modulation  type  classifier  rejects  the  input, 
the  EXEC  will  recover  provided  that  the  combination  classifier  and 


the  remaining  baud  or  type  classifier  agree 


The  Modem  ID  classifier  is  not  used  in  this  decision,  rather  it  is 
required  in  Phase  3  in  lieu  of  consistent  database  collection 
parameters. 


C.3.3  Phase  3  -  Channel  Impairment  Characterization 

Phase  3  consists  of  two  sub-phases;  the  statistical  channel  characterizer 
for  the  estimation  of  channel  S/N  and  harmonic  distortion,  and  the  phase 
jitter  detector. 

Programs: 


EVLIMP  the  peak-to-average  ratio  and  the  peak/average  ratio  variance  of 

the  rectified  signal  are  calculated  and  used  to  index  into  a  tem¬ 
plate  file  identified  by  the  modem  ID  classifier  of  Phase  2.  The 
requirement  for  the  fingerprint  level  ID  of  modem  is  imposed  by 
the  unexplained  discrepancies  in  the  signal  formats  due  most 
likely  to  the  inconsistent  collection  parameters,  sample  rate, 
etc. 


The  result  of  the  template  search  is  the  classification 
result  for  the  presence  and  degree  of  HD  or  kGN  on  the  line  in  the 
form  of  the  Phase  3A  executive  summary. 

is  the  IIR  filter  demodulator.  The  generic  identification  indi¬ 
cates  a  a^sk  file  (for  example,  DPSK1200.PRM)  which  contains  the 
mix  frequency  and  proper  filters  for  the  demodulation  of  this  gen¬ 
eric  type. 


DEMODI 


Fop  the  three  cases  of  bits/ symbol  1,  2,  and  3,  the  FM  tran¬ 
sition  removal  tasks  are: 

REMBIT 1 
REMBIT 2 
REMBIT 3 

The  actual  ID  of  the  modem  is  used  here  to  indicate  the  proper 
thresholds  for  the  operation.  The  thresholds  are  stored  on  disk 
under  the  name  S".PRM"  where  'S'  is  the  symbol  of  the  modem  sig¬ 
nal.  The  input  extension  ".FM"  is  assumed  and  ".NFM"  (No  FM)  file 
is  produced  for  specified  input  waveform  file. 

RESAMP  Applies  a  200Hz  low-pass  Butterworth  2-pole  filter  to  smooth  the 
".NFM"  input  file.  The  signal  is  resampled  by  4  and  stored  back  on 
disk. 

PWRSPT  As  described  earlier,  the  power  spectrum  generator  operates  on  the 

specified  input,  the  resampled  no-FM  signal  in  this  case,  and  pro¬ 
duces  a  single  line  power  spectrum.  This  is  plotted  for  the  user. 

MODPEAKDT  is  the  modified  peaks  detector  which  operates  on  the  FM  spectrum. 
The  output  is  as  shown  below: 


line 

(1)  M 

25Hz 

(2) 

7.22dB 

(3)  C 

OH  2 

(4)  «2dB,  OH z 

(5)  HP 

(6)  M 

19Hz 

(7) 

6. 14dB 

where 

line  (1) 

Is  the  maximua  frequency  component  in  the  specified  band 

line  (2) 

is  the  power  of  the  line  (1)  component 

line  (3) 

is  the  center  frequency  of  the  component  or  is  0  if 

unique  spike  is  found 

no 

line  (4) 

Is  the  component  bandwidth  at  2dB  below  the  peak,  or 

is 

0  if  no  component  is  found 

line  (5)  Is  the  cause  for  rejection  of  the  spike  (If  It  Is 
rejected) : 

HP  a  multiple  peaks 

PO  a  peaks  overlapping  band  edge 


\* 


line  (6)  is  the  maximum  frequency  of  the  second  highest  spectral 
component  in  the  specified  band 

line  (7)  is  the  maximim  power  level  of  the  second  highest 
spectral  component  in  the  specified  band. 

In  addition,  an  estimate  of  the  noise  floor  for  the  0-200Hz  band  is  given 
in  the  leftmost  column.  This  serves  as  a  reference  when  examining  spectral 
power  levels  of  the  peaks  detector  output. 


WORD 


DESCRIPTION 


1 

2,3 

4 

5 

6 

7 

8 
9 

13 

14 

17,18 


No.  of  Channels 
Sample  Rate  (1*4) 

Start  Time  (Hours) 

Start  Time  (Minutes) 

Start  Time  (Seconds) 

End  Time  (Hours) 

End  Time  (Minutes) 

End  Time  (Seconds) 

Run  Number 

Data  Type  (0;  Integer  9 
32-bit  Floating  Point) 

Number  of  Data  Records  (1*4) 


All  Words  are  1*2  unless  otherwise  specified.  Only 
the  sample  rate  and  data  type  are  used  by  the  MSA 
software. 

Table  C-l  IDRS  Waveform  Header 


c-ll 


B21D70064.PR 


Rectified  Signal  Power  Spectrum 


Power  Spectrum  of  2-stage  Lattice  Filter  Residue  Rectified 


PHASE  2:  EXECUTIVE  SUMMARY 


DATA  FILE:  B21D70064. VEC 
NUMBER  OF  DATA  VECTORS:  8 


LOGIC  FILE:  M0DTYP.F3K 
F  P  D  Q  * 

0  0  8  0  0 

MODULATION:  PSK/DPSK 


LOGIC  FILE:  BAUDRT . F3K 

1  2  3  4  5  • 

8  0  0  0  0  0 

BAUD  RATE:  1200 


LOGIC  FILE:  MODBRT . F3K 

FDdQqAP* 

08000000 

MODULATION,  BAUD  RATE:  PSK/DPSK, 1200 


LOGIC  FILE:  MODEM ID.F3K 

ebedhlaprei 
070000000  10 

MODEM  TYPE:  CODEX  LSI -48,  MODE  B 


FINAL  DECISION: 

MODULATION,  BAUD  RATE:  PSX/DPSI, 1200 


Figure  C-9 

Phase  2  Executive  Summary 


C-20 


J 


PHASE  3A  EXECUTIVE  SUMMARY 


REFERENCE  LIBRARY:  b 


PROCESSING  B21D70064.NRM 

•• -  RESULTS  OF  THE  PEAK/ AVERAGE  RATIO  CLASSIFIER  - •• 

••  NO  SIGNIFICANT  HARMONIC  DISTORTION  OR  GAUSSIAN  NOISE  DETECTED  •• 

•• - RESULTS  OF  THE  PEAK/ AVERAGE  RATIO  VARIANCE  CLASSIFIER - •« 

••  NO  SIGNIFICANT  HARMONIC  DISTORTION  OR  GAUSSIAN  NOISE  DETECTED  •• 
••••••••••  pHASE  3B;  EXECUTIVE  SUMMARY  •••••••••• 


FILE:  B21D70064 


NOISE 
4.05  db 


18.  hz 

60.  hz 

120.  hz 

180.  hz 

~H  23.  Kz 

nr- 

63.  hz 

r~T~ 

135.  hi  1  M 

i65.  hz 

7.22  db 

13.33  db 

5.04  db  ; 

l 

5.38  db 

C  0.  hz 

c 

60.  hz 

C 

0.  hz  !  C 

0.  hz 

9  2.dbi  0.hz 

9  2. 

db,  6.hz 

9  2. 

db,  0.hzl0  2. 

i 

db,  0.hz 

MP 

MP 

:  mp 

M  19.  hz 

M 

0.  hz 

M 

123.  hz  :  M 

195.  hz 

6. 14  db 

0.00  db 

5.04  db  : 

3.95  db 

Figure  C-10 

Phase  3  Executive  Summary 


C-21 


FM  Spectrum  After  Enhancement  Showing  60  Hz  Phase 
Jitter  at  64%  P/P  Enhanced  8dB  Above  the  Noise  Floor 


END 

DATE 

FILMED 

6—83 
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